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7 —). LAGINHA SERAFIM — Desastres de barragens em 1976. 
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O ensino da Geologia no 1.5. T. 
para o curso de Engenharia Civil. 
Evolução e visão prospectiva [ 


RESUMO 


Os autores apresentam um resumo da evolução 
do ensino da Geologia para estudantes do Curso de 
Engenharia Civil do I.8.T. nos últimos 20 anos. 

Em seguida é apresentada uma visão prospectiva 
do que deverá ser esse ensino e do equilibrio a haver 
entre os conhecimentos dos fundamentos da Geologia, 
das suas metodologias e das suas aplicações a pro- 
blemas concretos da Engenharia Civil, 

Às principais conclusões são apresentadas sob a 
forma de um diagrama onde se doseia a influência e 
à sobreposição que há entre os conhecimentos impor- 
tantes nas diferentes fases da resolução de problemas 
reais da Engenharia Civil fortemente dependentes da 
Geologia. 

São feitas, ainda, algumas considerações sobre os 
tipos de exames que os alunos realizam, 


| — INTRODUÇÃO 


Motivos vários têm contribuido para a incom- 
preensão do papel importante das Geociências no con- 
texto do Curso de Engenharia Civil. Nomeadamente, 
não é estranha a esse fenómeno certa má aceitação 
dos estudos de natureza geológica gerada, no liceu, 
por um ensino que hipertrofia a memorização. 

Acresce que a visão nebulosa que os alunos têm 
do interesse futuro das Geociências na sua actividade 


L, AIRES-BARROS e J. F. PINHARANDA 
Lab. de Minerologia e Petrologia, |. 5. T. 


ABSTRACT 


The authors present aq summary of the evolution 
of the teaching of Geology to Civil Engineering stu- 
dents at Instituto Superior Técnico (I.8.T.) during 
the last 20 years. 

After this they present their prospective vision 
on the teaching of Geology and Engineering Geology 
for that students. The aims of this teaching and the 
problems regarding the need of a balanced knowledge 
of the fundamentals of Geology, their methodologies 
and their applications to concrete problems of civil 
engineering are pointed out. The main bases conside- 
red in order to program the Geology course to Civil 
Engmeering students are presented in diagramatic 
form where scientific and technical fields required 
for the resolution of problem of engineering geology 
are refered., 

Some comments on examinations are also refered, 


profissional é outro óbice à sua boa aceitação na al- 
tura em que aparecem inseridas no curriculo escolar. 

Decorre do exposto que é prioritário destruir o 
preconceito anti-geológico com que os alunos se apre- 
sentam no início dos cursos sobre Geologia. Este tra- 
balho tem de ser acompanhado pelo despertar do 
interesse para uma ciência de âmbito diferente das 
restantes constituintes do curso de Engenharia Civil, 


(*) Comunicação apresentada na «HReunião sobre Investigação e Ensino em Geotecnia», realizada no L. N. E. €C. em 


18 de Dezembro de 1975. 
Manuscrito recebido para publicação em 6G/1/77. 
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ou seja, com características epistemológicas bem di- 
ferenciadas que exigem metodologia adequada. 

Conseguido este objectivo preliminar e fundamen- 
tal, hã que procurar transmitir conhecimentos basi- 
lares das Geociências, chamando a atenção para as 
suas implicações na concepção e execução das cbras 
de engenharia civil. 

Como fecho deste trabalho escolar e de dinami- 
zação e reformulação de conceitos, há que facultar 
aos alunos treino em algumas técnicas de uso habi- 
tual em Geologia de Engenharia. Deverão ser focados 
especialmente problemas ligados à interpretação cor- 
recta dos resultados, quantificados ou não. Deve ainda 
motivar-se o aluno para a escolha criteriosa das téc- 
nicas geológicas a usar e para a aceitação das suas 
eventuais limitações. 

Não tendo de modo algum a pretensão de formar 
especialistas neste domínio, tem sido preocupação 
constante promover a compreensão da terminologia 
geológica, indispensável ao entendimento e colabora- 
ção entre profissionais de engenharia civil e de geolo- 
gia aplicada. 

Este aspecto é de grande relevância como factor 
de êxito na realização e condução de obras em que 
estes dois domínios de conhecimento se aplicam. A 
experiência mostra, com clareza, que esta incompreen- 
são semântica está na origem de alguns inêxitos cuja 
responsabilidade se deve repartir igualmente pelos 
intervenientes, 


I— O ENSINO DA GEOLOGIA PARA O CURSO DE 
ENG." CIVIL ATÉ 1970 


No currículo escolar do curso de Engenharia Ci- 
vil que vigorou até 1955, figurava, no 3.º ano uma 
cadeira anual denominada «Curso Geral de Minera- 
logia e Geologia». Como disciplina facultativa, podiam 
os alunos frequentar «Geologia», cadeira do 4.º ano 
do Curso de Engenharia de Minas. 

A reformulação curricular operada em 1955 nada 
de novo trouxe quanto ao ensino das cadeiras refe- 
ridas para o Curso de Engenharia Civil. Para sermos 
mais precisos, mudou o nome do «Curso Geral de Mi- 
neralogia e Geologia» para «Mineralogia e Geologia 
Gerais»! Esta mudança de nome não foi tão inofensiva 
como pode parecer... É que o diploma de 1955 a que 
nos referimos (Decreto n,º 40 378 de 14 de Novembro 
de 1955) criou uma cadeira com o mesmo nome, para 
o Curso de Engenharia Quimico-Industrial (actual 
curso de Engenharia Química), situada no 1.º ano. 
Este facto levou algumas Escolas Superiores a fundir 
os dois cursos (1.º ano para «Química» e 3.º ano para 
«Cívil ), dando um programa comum! 

No I.S.T. conseguiu-se evitar esta fusão que a 
burocracia aconselhava, mas o bom senso repudiava. 
Todavia não foi pacífica a tentativa (com êxito) do 
desdobramento de «uma mesma disciplina», em duas. 

Os princípios que informaram esta cadeira («Mi- 
neralogia e Geologia Gerais») podem considerar-se 
bastante «tradicionais» e um pouco «enciclopedistas», 

Diremos «tradicionais» na medida em que se 


procurou ajustar a uma cadeira anual e «geral» o 
essencial de Geologia Geral, seguindo a experiência 
pessoal bastante apoiada na literatura, aliás abun- 
dante, Diremos ainda um pouco «enciclopedistas» por- 
que, com critério e conhecendo as dificuldades, se 
optou por iniciar o curso pelos princípios gerais da 
Cristalografia, passando-se aos elementos da sistemá- 
tica mineralógica com ênfase para os silicatos, Sobre 
esta base se apoiavam as lições de petrografia. Nisto 
se consumia metade do curso. A seguir expunham-se 
os fundamentos da Geologia Geral (Geomorfologia, 
Tectônica, Elementos de Estratigrafia) e um pouco 
de Geologia de Portugal, Dava-se certa ênfase a ca- 
pítulos como a estrutura dos silicatos, os minerais 
argilosos, os fenómenos de alteração e alterabilidade 
de rochas, Faziam-se algumas excursões aos arredo- 
res de Lisboa (particularmente à região de Cascais 
e Sintra) e nos últimos 3 anos de regência desta ca- 
deira fez-se uma excursão mais longa, de vários dias, 
cobrindo o vale do Tejo, o percurso de V. Velha de 
Ródão até à Serra da Estrela, Em seguida, das re- 
giões do Maciço Hercínico descla-se pelo Buçaco até 
à Orla Mesozoica ocidental e, por Coimbra e Porto de 
Mós, penetrava-se no maciço calcário estremenho para 
regressar a Lisboa pela margem direita do Tejo. 

Os cursos práticos acompanhavam «pari pasasu> 
o curso teórico. Dava-se particular ênfase à petro- 
grafia (análise macroscópica) e ao estudo das cartas 
geológicas. A culminar, havia que realizar um corte 
geológico a partir da carta, interpretá-lo com certo 
pormenor e elaborar uma memória descritiva. 


A traço grosso foi este o «ambientes de uma 
cadeira frequentada anualmente por cerca de mela 
centena de alunos, Cadeira aceite, na generalidade, 
com agrado pelos alunos, até pela feição prática, com 
chamamento à observação sistemática, com excursões 
ao terreno e culminando com uma excursão de alguns 
dias ao centro do país. Esta excursão foi sempre or- 
ganizada integralmente pelos alunos embora coadju- 
vados em tudo pelo corpo docente. 


Quanto à cadeira facultativa de Geologia, era 
frequentada por uma percentagem diminuta de alu- 
nos, 2 a 5% (!), sempre os melhores estudantes, AÍ 
se lhes dava uma feição mais aplicada aos princípios 
que levavam do curso geral, 


HWi— O ENSINO DA GEOLOGIA PARA O CURSO DE 
ENG.* CIVIL A PARTIR DE 1970 


O Decreto n,º 640/70 de 10 de Novembro de 1970 
procurou proceder à actualização dos planos de estu- 
dos dos curzos de engenharia professados nas Univer- 
sidades portuguesas «de harmonia com as normas 
pedagógicas que a experiência tem revelado eficien- 
tes» como se afirma no seu preâmbulo, 

Para os cursos professados no I.S.'T. no domí- 
nio do curso de Engenharia Civil, optou-se por um 
curso de «Mineralogia e Geologia Gerais» no 1.º se- 
mestre do 2.º ano e por um curso de «Geologia Apll- 
cada» no 2.º semestre do 3.º ano, 
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Desde já se fazem alguns reparos a este posi- 
cionamento das mencionadas disciplinas. 

a) Tem-se revelado contraproducente a separação 
das duas cadeiras por um semestre em «branco», 


b) É fundamental que a Topografia preceda o 
curso de «Mineralogia e Geologia Gerais». 


c) É controverso o início deste curso eminente- 
mente prático no 2.º ano ao lado de disciplinas de 
indole muito diferente, como a Física e a Matemá- 
tica, e, sobretudo, quando os alunos ainda não dis- 
põem de uma perspectiva da sua futura actividade 
profissional, 


Do exposto decorre que estas duas cadeiras, a 
manterem-se, numa reformulação futura deviam co- 
locar-se: a de «Mineralogia e Geologia Gerais» no 
1.º semestre do 3.º ano (trocando com Topografia); 
a de «Geologia Aplicada» deve permanecer no 2.º se- 
mestre do 3.º ano. 

Vejamos agora como foi encarada a reformulação 
do curso anual referido no capítulo anterior à luz dos 
novos interesses, tendo em atenção o aparecimento 
de um semestre dedicado à «Geologia Aplicada», 

No primeiro ano de aplicação do novo esquema 
curricular foi esboçada a tentativa de: 

a) condensar na cadeira semestral de «Minera- 
logia e Geologia Gerais» o essencial do curso anual 
anterior, 

b) reservar a nova cadeira semestral de «Geolo- 
gia Aplicada» para leccionar matérias mais directa- 
mente ligadas a problemas de geologia de engenharia, 

Logo se verificou que tal via só se poderia con- 
cretizar à custa dum acentuado abaixamento do nível 
e pormenor do curso, 

De facto, sem tal «sacrifício», haveria «trans- 
bordo» de matéria da cadeira geral para o curso de 
Geologia Aplicada, o que, nas circunstâncias então 
existentes, foi considerado preferível, 

Mas compreendeu-se a necessidade de uma re- 
formulação mais radical do curso, de molde a per- 
mitir ajustar as duas cadeiras semestrais a matérias 
mais conzentâneas com as respectivas designações, 
isto é, reservando de facto a cadeira de «Geologia 
Aplicada» para questões de geologia de engenharia. 
Essa profunda modificação do curso não poderia ser 
pfectivada de imediato. 

Assim, nesses primeiros tempos, a cadeira se- 
mestral de «Mineralogia e Geologia Gerais» foi preen- 
chida, segundo o esquema clássico, por noções de 
cristalografia, mineralogia (com ênfase na mineralo- 
gia das rochas e breves referências aos minerais me- 
tálicos) e petrografia, incluindo alteração e alterabi- 
lidade das rochas. 

Na cadeira semestral de «Geologia Aplicada», a 
matéria abrangia assuntos relacionados com Geologia 
Fisica e um breve bosquejo da Geologia de Portugal. 
Como conhecimentos de Geologia de engenharia, uma 
introdução aos métodos de prospecção, com especial 
ênfase para o estudo das sondagens, algumas lições 
de Hidrogeologia, e algumas técnicas de ensaios geo- 
técnicos. 

Nas aulas práticas, os alunos acompanhavam, 
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nos moldes habituais, a matéria leccionada nas aulas 
teóricas, 

Vicissitudes várias, nomeadamente interrupção 
das aulas por movimentos reinvidicativos, impediram 
que até ao ano lectivo de 1974/75, se concretizasse 
o objectivo acima referido de concentrar toda a ma- 
téria introdutória na cadeira de «Mincralogia e Geo- 
logia Gerais», permitindo alargar o campo de ensino 
da Geologia de engenharia. 


IV — LINHAS GERAIS PARA A ORIENTAÇÃO DUM PRO- 
GRAMA DAS CADEIRAS DE GEOLOGIA PARA O 
CURSO DE ENG. CIVIL 


Referiu-se que os objectivos da inclusão das ca- 
deiras de âmbito geológico mo curso de Engenharia 
Civil não poderiam ser os de formar técnicos espe- 
cializados em Geologia Aplicada, mas, primordial- 
mente: 


— chamar a atenção para a necessidade de enca- 
rar as implicações geológicas na concepção e 
execução das obras de engenharia civil; 

— obter um mínimo de compreensão da termi- 
nologia geológica, para facilitar o entendimen- 
to e a colaboração entre o engenheiro clvil c o 
especialista de geologia aplicada. 


Esta limitação de objectivos não resulta da ccn- 
sideração de incompatibilidade entre o engenheiro ci- 
vil e qualquer especialização no domínio da geologia 
aplicada, mas somente porque essa eventual via de 
actuação não se poderá considerar normal como «sai- 
da» dum curso de engenharia civil, 

Para atingir aqueles objectivos parece acertada 
a realização de alguns trabalhos práticos, em parte 
de carácter mais acentuadamente geológico, a partir 
dos quais se procurarão focar aspectos de aplicação 
da geologia à engenharia civil, nomeadamente por 
discussão de casos concretos e por apreciação de tra- 
balhos realizados, e também pelo estudo de bases 
técnico-científicas indispensáveis para a resolução das 
questões e solução de problemas de aplicação. 

O problema das matérias que deverão ser foca- 
das no decurso da cadeira de «Geologia Geral» não 
pode ser definido sem ser em conjugação com o da 
cadeira complementar de «Geologia Aplicada», e sem 
considerar os cursos de «Mecânica dos Solos» e de 
«Mecânica das Rochas», que fazem parte — obrigató- 
ria ou facultativamente — do curso de Engenharia 
Civil. 

Para o estudo desze programa, parece possível 
classificar os conhecimentos de índole geológica no 
que respeita à sua aplicação em engenharia civil, em 
três domínios interdependentes e interpenetrados: 


— Bases cientificas fundamentais 
— Métodos gerais de estudo e ensaio 
— Aplicações a obras de engenharia civil 


Com vista a delinear um programa de estudos 
adequados aos interesses do futuro engenheiro civil 


e do geólogo de engenharia no decurso da resolução 
de problemas genéricos. Isto será feito com ajuda do 
esquema apresentado adiante, que se procurará ex- 
plicar. (Fig. 1). 


A primeira etapa será a definição do problema, 
da responsabilidade dos técnicos de engenharia civil, 
para a qual se deverão servir dos seus conhecimentos 
dos problemas específicos de geologia aplicada ao 
género de obra em estudo e das bases científicas no 
domínio geológico. 

A partir dessa definição, os técnicos de geologia 
de engenharia, em colaboração com os de engenharia 
civil, nomeadamente os especialistas de mecânica das 
rochas e dos solos, elaborarão o programa de reco- 
nhecimento. Para este são essenciais profundos co- 
nhecimentos de processos de prospecção e de ensaio 
de geologia aplicada ao tipo de obra e das bases 
científicas fundamentais. A colaboração dos técnicos 
de engenharia civil é essencial para a definição de 
prioridades na execução do programa, de forma a que 
os primeiros elementos a obter possam dar indicações 
decisivas para a localização e concepção da obra e 
possam ser aproveitados para eventuais ajustamentos 
do resto do programa de prospecção. 

A campanha de reconhecimento, em geral a cargo 
de uma empresa especializada, e estreitamente acom- 
panhada pelos especialistas de geologia de engenha- 
ria, obriga, principalmente, a boa compreensão dos 
processos de reconhecimento. Se bem que neste esca- 
lão, pareça secundário o papel dos técnicos de enge- 
nharia civil, considera-se de grande utilidade que 
estes acompanhem com atenção o desenrolar da cam- 
panha, 

Outra fase que se segue obriga a nova colabora- 
ção entre geólogos de engenharia e engenheiros civis: 
trata-se da análise dos resultados, para a qual é fun- 
damental um profundo conhecimento das caracteris- 
ticas específicas dos vários métodos de investigação 
utilizados, das suas virtualidades e limitações no que 
respeita ao caso concreto em estudo, das bases cienti- 
ficas para estruturar interpretações correctas diver- 
sificadas e para as discutir, e ainda das questões de 
geologia aplicada às obras. Esta fase deve acompa- 
nhar a própria execução da campanha de reconheci- 
mento. 

O escalão seguinte é o da apresentação de um 
relatório de reconhecimento com as indicações, infor- 
mações e recomendações para a concepção e execução 
das obras. Muito frequentemente, tanto o geólogo de 
engenharia como o engenheiro civil, admitem que 
este relatório fecha a intervenção do primeiro. Tal 
ideia parece errada pois, por melhor que seja o espe- 
cialista de geologia aplicada e por mais fiéis que 
tenham sido os resultados da campanha de reconhe- 
cimento, a sua informação final resulta sempre duma 
interpretação, duma generalização a partir de ele- 
mentos dispersos, pontuais ou lineares, quanto muito, 
para a caracterização do volume interessado pelas 
obras. Portanto, o geólogo de engenharia tem de es- 
clarecer as suas informações na fase de elaboração 
do projecto, na qual os técnicos de engenharia civil 


aplicarão os seus conhecimentos de geologia aplicada 
e de geologia geral discutindo com aquele as impli- 
cações que a sua interpretação acarreta na concepção 
do trabalho. Finalmente, o especialista de geotecnia 
deverá acompanhar a execução da obra, em inteira 
colaboração com os engenheiros do estaleiro, para 
ajustar a sua interpretação à medida que forem sur- 
gindo elementos fornecidos pelos próprios trabalhos 
e para orientar as necessárias operações de consoli- 
dação, e tratamento dos terrenos e mesmo de reco- 
nhecimento complementar, 

Em resumo, os técnicos de engenharia civil, no 
que diz respeito a trabalhos relacionados com geolo- 
gia aplicada, intervém principalmente: 

— na definição do problema, tendo em atenção as 
características das obras e suas implicações 
epológicas; 

— na programação dos trabalhos de reconheci- 
mento, definindo prioridades e ensaios geotéc- 
nicos necessários; 

— ma apreciação dos resultados e na adaptação 
das obras às características reveladas pelo re- 
conhecimento. 


Embora estas intervenções sejam em colaboração 
com geólogos de engenharia, elas obrigam os enge- 
nheiros civis a utilizar técnicas e conhecimentos geo- 
lógicos e geotécnicos. 

Assim, no decurso das duas cadeiras de âmbito 
geológico do curriculo de Engenharia Civil parece 
indispensável focar diversos capítulos dos três domi- 
nios de conhecimento considerados, com especial aten- 
ção para os métodos de prospecção e para as bases 
científicas — que se pensa poderem ser preferencial- 
mente objecto de estudo teórico-prático — referindo 
sempre que possível, nas mesmas aulas teórico-práti- 
cas e também em aulas de tipo conferência aberta à 
discussão, como é óbvio, aspectos ligados aos pro- 
blemas especiais de aplicação. Sempre que possível, 
estes conhecimentos teóricos ou teórico-práticos, de- 
verão ser reforçados em visitas de estudo. 

É impossível leccionar um curso completo de 
Geologia Aplicada à Engenharia Civil em duas ca- 
deiras semestrais, e pensa-se mesmo que não é esse 
o objectivo da respectiva inclusão no currículo lectivo. 
Hã necessidade, portanto, de optar pelo maior desen- 
volvimento de determinados assuntos em prejuízo dos 
restantes, 

Nos campos das bases científicas e dos métodos 
de reconhecimento e ensaio parece, contudo, preferi- 
vel referir nas aulas uma gama bastante geral dos 
assuntos, embora sem grande profundidade, incompa- 
tível aliás com a especialização dos alunos. Esta re- 
lativa «leveza» não pode significar a redução do cur- 
so a um somatório de definições, pois julga-se impres- 
cindível a verdadeira compreensão dos fenómenos e 
das técnicas. 

Já no domínio da Geologia Aplicada à Engenha- 
ria Civil, se pensa vantajoso estabelecer, de acordo 
com os alunos, um esquema de prioridades para de- 
senvolvimento dos assuntos a leccionar, necessaria- 
mente em número limitado, 
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V— O PROBLEMA DE AVALIAÇÃO DOS CONHECI- 
MENTOS 


Até 1970, o sistema de avaliação de conhecimen- 
tos seguiu o método tradicional do exame escrito e 
oral individual, O exame oral era requerido aos alunos 
cuja informação da prova escrita se situava no do- 
mínio 9-12 valores. Abaixo daquela informação o 
aluno estava reprovado, a menos que a informação 
prática fosse relevante. Neste caso o aluno era admi- 
tido a exame oral alargado e profundo. Acima de 12 
valores o exame oral era facultativo, 

Todavia nos anos lectivos 1969/70 e 1970/71 en- 
salou-se, a solicitação dos alunos, o exame colectivo. 
Este exame versava sobre temas propostos e desen- 
volvidos por grupos de alunos (de 3 a 5) em trabalho 
escrito. Estes documentos, bem como os trabalhos 
práticos, em especial o corte geológico e seu relatório, 
eram analisados em exame oral colectivo. O júri era 
constituído por 3 elementos que trocavam «impressões 
examinatórias» com o grupo de discentes, 

Dois problemas se levantaram: 


1 -—A despeito de se apresentar um vasto leque 
de assuntos para desenvolvimento, houve polarização 
do interesse dos alunos em um núcleo restrito de 
temas. Daqui resultou certa repetição de trabalhos. 


2 — Com o empolamento do número de alunos 
tornou-se difícil em tempo oportuno «ver e corrigir» 
o grande número de relatórios. Tenha-se em atenção 
que o espacejamento inter-exames das várias disci- 
plinas era de 8 a 10 dias. Daqui a dificuldade prática 
de corrigir e examinar oralmente tão largo número 
de grupos de trabalho. 


Como processo de avaliação de conhecimentos ve- 
rificou-se que era demasiado elevado o número de 
alunos que «vivia à sombra» dos 2 a 3 trabalhadores 
do grupo. Deste modo, cerca de 50% dos componentes 
dos grupos era transportado por arrastamento. 

Todas as experiências quer de transmissão de 
conhecimentos, quer da sua avaliação tem deparado, 
nas cadeiras em causa, com a questão dos grandes 
números. Os signatários têm experiência colhida quer 
em cursos de pequenos números (v. g. Minas e Me- 
talurgia), quer dos grandes números de alunos (v. E. 
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Civil e Química) e ainda da expansão do curso de civil 
(passando de 50 a 200 alunos) com não acompanha- 
mento de dilatação correlativa da estrutura (instala- 
ções, material, docentes, etc.). 


Com a reformulação das referidas cadeiras, tam- 
bém houve que repensar o sistema de avaliação dos 
conhecimentos. Assim adoptou-se quer o sistema dos 
testes sintéticos com pergunta e respostas múltiplas 
competindo ao aluno eleger a afirmação certa, quer 
os testes teórico-práticos com relativa frequência ao 
longo do período escolar. Estas experiências foram 
positivas, apenas havendo a referir que há que baixar 
drasticamente o número de alunos por docente a fim 
de não tornar estes máquinas de testar alunos. Refi- 
ra-se ainda a dificuldade de organizar testes signifi- 
cativos, onde os alunos possam, de facto revelar os 
seus conhecimentos, 


O ideal seria quase abolir o teste formal. Para 
atingir uma forma racional de avaliação de conhe- 
cimentos há que aliar a razoável razão discente/ 
(docente com uma melhor consciência pedagógica 
dy aluno médio português. Esta, aliás, decerto se 
fortaleceria com a melhoria da relação, docente/ 
“discente, que talvez possibilitasse formas pedagogi- 
camente correctas de isolar e recuperar — ou rejeitar, 
em casos extremos — o cábula, o oportunista, o para- 
sita, dos que estudam. E este que, defraudando de 
todas as formas o país, procurando caminhar na 
sombra dos mais honestos ou mais válidos, se agarra, 
quando o não consegue, a um tipo degradante de 
contestação, que mais não é que a defesa cega, e a 
todo o custo, de privilégios sociais injustos. 


Cabe aqui afirmar o respeito que merece, como 
sempre mereceu nestas cadeiras, a contestação vivi- 
ficadora, que não procura «facilidades», mas a mais 
correcta orientação do esforço colectivo docentes /dis- 
centes que deve ser o trabalho do ensino, 


Mais ultimamente a avaliação dos conhecimentos 
diluiu-se e as suas consequências devem ser conside- 
radas. Todavia ensaiou-se a discussão oral dos tra- 
balhos executados nas aulas práticas pelos alunos, 
agora em grupos de 2 e excepcionalmente de 3, ca- 
bendo aos examinandos a «defesa», isto é, a justifi- 
cação das soluções apresentadas, 
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RESUMO 


Este trabalho trata da rotura de barragens 
durante o ano de 1976 e está na sequência do tra- 
balho apresentado pelo autor sobre o mesmo assunto 
no Laboratório Nacional de Engenharia Civil em 
1960. Apresentam-se descrições de roturas em bar- 
ragens relativamente pequenas (mais de lôm de 
altura) na Carolina do Norte, no Paquistão e na 
Noruega, mas o trabalho trata essencialmente da 
rotura da barragem de Teton nos Estados Unidos, 
em Junho de 19%. Esta é a barragem mais alta 
(1083m) que até hoje rompeu. Apresenta-se a sequên- 
cia dos acontecimentos durante a rotura, atribuindo 
esta ao facto da rocha ser muito permedvel, à inefi- 
ciência das cortinas de injecção, e ainda ao facto de 
existir uma trincheira corta-águas, na qual se pro- 
duziu um efeito arco que possivelmente permitiu 
fracturação hidráulica. A erosão interna e subse- 
quente erosão no aterro perto da rocha foram, sem 
dúvida, as causas principais. 


1. ROTURAS EM 1976 


O presente trabalho diz respeito às roturas de 
barragens que tiveram lugar durante o ano de 1976. 

Foi em 1960 que, também no Laboratório Nacio- 
nal de Engenharia Civil, tivemos ocasião de referir 
os acidentes graves em barragens ocorridos em 1959. 
Terão, todos aqueles que estão interessados nessas 
estruturas, tido ocasião de notar que esse também 
foi um ano tristemente célebre na engenharia. De 
facto, dão-se então acidentes, únicos na história da 
construção de barragens, em França, em Espanha 
e no Brasil, onde, efectivamente, romperam obras 
importantes, uma de betão, outra de alvenaria e outra 
d: terra e também se produziram alguns graves aci- 
dentes que não chegaram até à rotura, mas que 
trouxeram igualmente graves preocupações, na Índia, 
no Peru e na Jugoslávia. Portugal, poucos anos de- 
pois, não escapou à «lista negra», vendo a destruição 
por uma cheia de uma barragem de enrocamento 
quase concluida. 


A rotura de uma barragem, principalmente quan- 
do ela é importante, é um acontecimento de tal forma 
grave e insólito que, quase sempre, muda o caminho 
das ideias e os critérios de projecto destas obras, 
De facto, sempre tem sido assim. Neste ano, mais do 
que nunca, é de crer que tenhamos que rever muitas 
das ideias expressas nos tratados. 1976 foi um ano 
tristemente célebre, na medida em que a imprensa 
diária e a imprensa técnica noticiaram seis acidentes. 
Porém, quanto aos casos que óocuparam os jornais, 
o da Colômbia não se passou com uma barragem, 
mas com uma outra obra hidráulica. Foi um canal 
de encosta que rompeu pelo lado exterior, tendo a 
água enchido uma depressão que havia nessa encosta, 
o que produziu repentinamente o seu desabamento, 
indo a água armazenada e as rochas destruir grande 


SUMMARY 


This paper concerns the rupture of dams during 
1976 and continues the work presented by the author 
in 1960 at the LNEC (National Laboratory of Cívil 
Engineering) on the same subject. Descriptions of 
ruptures of relatively small dams (above 15m) im 
North Carolina, Pakistan and Norway are presented. 
Most of the paper concerns the rupture of the Teton 
Dam in the United States in June, 1976. This was 
the highest dam (108m) to ever have failed. The 
sequence of events is presented with photographs. 
The rupture is attributed to the very pervious rock, 
the inefficiency of the grout curtains and the exis- 
tence of a key trench, in which arching was produced, 
possibly allowing hydraulic fracturing. Piping and 
subsequent intense erosion of the soil near the rock 
were certanly major causes, 


parte da aldeia de Pereira, na Colômbia, onde, pelo 
menos, 80 pessoas desapareceram do mundo dos vi- 
vos. A Colômbia tem sido um país muito castigado 
com estes acidentes. Já em 1972 tinha caido lá uma 
outra obra que provocou um número de mortos se- 
melhante. Também os jornais e revistas técnicas 
anunciaram a queda de uma barragem no Vietnam. 
Hã, no entanto, lugar para dúvidas sobre se teria 
efectivamente sido uma rotura. Na Carolina do Norte 
(E. U. A.), no Paquistão e na Noruega houve efecti- 
vamente rotura de barragens que, tal como a que nos 
vai ocupar mais, a barragem de Teton nos Estados 
Unidos, são todas de terra e com mais de quinze 
metros de altura, 


A barragem de Carolina do Norte, com 17 m de 
altura, era uma obra de uma entidade particular, 
pelos vistos construída com pouca técnica. Em 1963 
já se notava que todo o seu paramento de jusante 
aparecia saturado, pelo que estava sujeita a certa 
vigilância por parte das autoridades desse Estado. 
Mesmo assim, não foram a tempo com o remédio e a 
barragem caíu por si própria, depois de chuvas in- 
tensas, tendo provocado quatro mortes. Pelos elemen- 
tos que dispomos, talvez se possa chegar à conclusão 
que houve erosão interna, 

A barragem de Bolan, no Paquistão, com 19m 
de altura, 50m de comprimento, 800 000 m3 de terras 
e 99 000 000 m3 de água armazenada, foi completada 
em 1961 e foi totalmente destruída por uma grande 
cheia que ocorreu em Setembro. Essa cheia arrastou 
e rompeu, provavelmente, mais oito pequenas barra- 
gens (todas de altura inferior a 15m portanto). 
Os efeitos desta catástrofe, onde avultam principal- 
mente os devidos ao carácter excepcional da cheia, 
que a rotura das barragens ampliou, teriam segura- 
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mente sido menores se o descarregador da barragem 
de Bolan tivesse maior capacidade. De facto, as noti- 
cias indicam que o descarregador foi absolutamente 
insuficiente para todo o caudal e a água passou sobre 
a barragem, destruindo-a completamente. Foram 
muitas as aldeias afectadas, donde 10 000 habitantes 
foram evacuados por helicópteros e outros meios, 
tendo desaparecido algumas centenas de pessoas. 


Na Noruega caíu em Março a barragem de 
Rappojern, uma obra hidroeléctrica projectada com 
boa técnica por uma firma consultura que já havia 
usado o mesmo tipo de projecto em dez outras obras 
semelhantes. Não conseguimos até hoje obter nenhu- 
mas informações mais. O assunto está em inquérito 
na Ncruega e estão prometidas, a todo o mundo 
teécnico, explicações pormenorizadas. O que não há 
é dúvida que foram vários os factores intervenientes 
na rotura desta barragem da Noruega, entre eles a 
erosão interna ao longo da descarga de fundo que 
foi montada numa estrutura em betão deixada no 
corpo do aterro. 


Se somarmos a estes casos atrás mencionados, 
o caso da barragem de Tarbela, no Paquistão tam- 
bém, que, desde 1974 estã em série risco devido a 
sucessivos acidentes nos túneis das descargas de fun- 
do e nas suas comportas, nas respectivas bacias de 
diss'pação e ainda na própria barragem, acidentes 
esses que têm motivado que a albufeira da obra não 
tenha sido totalmente cheia até 1976, temos que con- 
cluir que existem efectivamente muitas preocupações 
no mundo inteiro acerca deste assunto. Assim é que 
a Fundação de Engenharia dos E. U. A. («Engineering 
Foundation») e sociedades americanas de engenheiros 
promoveram em Dezembro, em Asilomar, na Cali- 
fórnia, uma reunião sobre segurança de barragens. 


2. SISTEMAS DE OBSERVAÇÃO DE BARRAGENS 


Em face dos desastres ocorridos em 1976, é oca- 
sião para dizer que, até hoje, nenhum sistema de 
auscultação das barragens que romperam deu notória 
indicação de perigo ou o alarme de que o acidente 
se ia produzir, Em geral, esses sistemas são de lei- 
tura manual e as observações são efectuadas com 
grande intervalo (semanas ou até meses), sendo, na 
maior parte das vezes, os cálculos e os diagramas de 
interpretação feitos manualmente. Perante esse facto, 
um dos pontos que tem sido muito discutido tem sido 
o interesse, e a necessidade até, de colocar sistemas 
automáticos de observação e controlo da segurança 
das barragens e de alarme às populações a jusante. 
Fizemos uma proposta, há dois anos, à nossa Comis- 
são Nacional de Barragens, proposta que foi até ao 
executivo da Comissão Internacional (CIGB), para a 
criação de um comité técnico que estudasse e propu- 
sesse à (CIGB) normas sobre esses sistemas em face 
da experiência actual. Infelizmente, não vimos ainda 
aprovada essa proposta, mas os engenheiros de gran- 
des barragens são unânimes em que se desenvolvam 
a automatização da auscultação e os sistemas de 
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alarme dessas obras. De facto, parece que aqui esta- 
mos atrasados relativamente a outros empreendimen- 
tos de grande porte: qualquer central térmica moder- 
na tem um sistema automático de observação e de 
alarme, tal como as centrais atômicas, evidentemen- 
EO.» 


3. A BARRAGEM DE TETON 


Na história da rotura de barragens, a obra mais 
alta que havia falhado era a de Baldwin Hills, tam- 
bém de terra, nos E. U. A., com 80m de altura, isto 
quando a barragem mais alta, já em construção, 
atingirá os 320 m. Neste contexto, a rotura da barra- 
gem de Teton, de que vamos agora tratar, é um 
dos acidentes mais importantes da história dos de- 
sastres na Engenharia Civil e é um caso absoluta- 
mente à parte. Além de ser a barragem mais alta que 
até hoje caíu (103 m), a rapidez com que se processou 
o acontecimento foi impressionante. Desde que a ro- 
tura era evidente até que ela se produziu, não se 
passaram sequer três horas. O acidente pôde ser 
fotografado e filmado por vários assistentes. O Bu- 
reau of Reclamation pôs à nossa disposição uma 
série de diapositivos e, depois de realizada esta con- 
ferência, o relatório da comissão de especialistas no- 
meada pelo Ministério do Interior do Governo Fede- 
ral Americano e pelo Governador do Estado de Idaho, 
que adiante referimos. O Bureau of Reclamation é a 
entidade americana governamental que até hoje cons- 
truiu maior número de barragens. Esse organismo, 
fundado no fim do século passado, tinha até hã bem 
pouco tempo os «records» americanos e mundiais da 
altura, tanto em barragens de betão como nas de 
terra. Projectou e construiu para cima de 260 barra- 
gens de terra e nunca teve an“*eriormente nenhum 
acidente com uma barragem sua. A barragem tinha 
um comprimento de 950 metros (fig. 1) e um volume 
d: oito milhões de metros cúbicos de materiais. Era, 
pelo que foi dito, uma barragem de terra, compactada 
mecanicamente, conforme perfil indicado na fig. 1, 
com um núcleo impermeável bastante espesso em 
argila e silte e espaldas de materiais aluvionares 
permeáveis, protegidas por enrocamentos. O reserva- 
tório teria uma capacidade de trezentos e cinquenta 
milhões de metros cúbicos e, quando a barragem 
rompeu, já tinha armazenados duzentos e cinquenta 
milhões de metros cúbicos. A barragem em si custou 
dois milhões de contos. Foi construída pela firma 
Morrison-Knudsen em associação com a Kiewit. 


As fundações da barragem são constituídas por 
uma rocha vulcânica (riolito ou tufo vulcânico) que, 
em certas zonas da margem esquerda, estava depo- 
sitado sobre material aluvionar recente, o que quer 
dizer que o derrame vulcânico foi muito recente, Ha- 
via muitas preocupações com esse aluviões subjacen- 
tes ao tufo vulcânico, mas o facto é que a rotura se 
deu na outra margem, onde tal situação não ocorre. 
Esta situação de lavas sobre aluviões quaternárias 
é relativamente comum, como já tivemos ocasião de 


verificar numa zona de vulcões na Costa Rica. No 
meio do vale havia uma zona de aluviões actuais 
que, em parte, foi escavada (fig. 2). O tufo vulcã- 
nico estava muito fracturado, cheio de diaclases mui- 
to abertas e portanto tinha uma permeabilidade mui- 
to alta. Só depois de ensaios de injecções de imper- 
meabilização é que foi decidido construir a barragem. 
Em 1973, quando a barragem já estava em início de 
execução, um grupo de geólogos do Geological Survey 
fez a afirmação (aliás suportada por elementos do 
célebre e reaccionário «Sierra Club» e outros grupos 
de conservacionistas) que a barragem não tinha con- 
dições de segurança, porque estava numa zona alta- 
mente sísmica. Isso fez com que fosse talvez, até 
hoje, das obras onde foram colocados maior número 
de sismógrafos. Bem analisados os argumentos que 
então se produziram, e estão todos escritos, chega-se 
à conclusão que aquilo que foi dito não tem nada 
que ver com o que aconteceu. Mas o facto é que há 
alarme mundial acerca dos sismos. Creio que esse 
alarme relativamente à acção sísmica não tem o 
fundamento que a experiência ou os cálculos deve- 
riam confirmar. Todavia o problema dos sismos não 
é de forma nenhuma um problema pouco importante. 
Direi é que pela maneira como, em matéria de es- 
forços de investigação, as coisas estão seguindo, se 
estão esquecendo problemas essenciais da segurança 
das barragens que merecem muitos esforços de pes- 
quiza. 

Na fig. 1 apresenta-se uma planta da barragem 
e respectivos órgãos acessórios, que incluem um des- 
carregador, na margem direita, uma pequena central 
com três grupos geradores na margem esquerda, uma 
tomada de água e também uma descarga de fundo 
que aproveitaram a galeria de desvio e ainda uma 
descarga de fundo auxiliar na margem direita por 
debaixo do descarregador. Na mesma figura apresen- 
ta-se um perfil transversal da barragem pela zona 
mais profunda do vale e também uma secção lateral, 
fazendo notar as trincheiras corta-águas, tanto a da 
rocha como a do aluvião (pormenores A e B). 

A fig. 2 mostra a planta e o perfil longitudinal 
das escavações onde se apreciam essas trincheiras, 
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tanto no aluvião na parte inferior do vale, como as 
laterais na rocha. Há que assinalar o perfil das trin- 
cheiras desses corta-águas escavadas nas encostas 
rochosas das duas margens. Essas trincheiras, com o 
perfil muito abrupto, têm aproximadamente 20m de 
profundidade e 12m de largura, Em baixo, e a partir 
do fundo das trincheiras, foram feitas três cortinas 
de injecção, as quais seguiram todos os métodos que 
é costume aplicar ao projectar e executar tais corti- 
nas, entre eles o da injecção de furos sucessivos por 
um processo «try and cut»: quando os furos recebem 
facilmente calda de injecção, perfura-se e injecta-se 
a uma distância intermédia uma nova linha de furos. 
Temos também a este respeito um pequeno pormenor 
que tem o seu interesse, e que é a protecção do início 
dos furos centrais com uma capa continua de betão 
(ver fig. 2). 

As duas condutas de descargas de fundo estão 
representadas nesta figura. Faz-se notar que, desde 
hã alguns anos, o critério do Bureau of Reclamation, 
e aliás de algumas outras entidades, é o de duplicar 
os órgãos de esvaziamento. Creio que este ponto ficou 
bem claro depois das discussões sobre acidentes em 
barragens por volta dos anos 62 e 63 e também de- 
pois da celeuma produzida pela primeira RANKINE 
LECTURE em Londres e das dificuldades por que 
passaram certas obras que não tinham descargas 
de fundo. Se Tarbela não tivesse descargas de fundo, 
e potentissimas, provavelmente já lá teria havido coi- 
sas bem piores do que as que já ocorreram. 

A fig. 3 mostra a obra já na fase final da cons- 
trução, pouco tempo antes do acidente. Ela ficava 
situada no rio TETON, um dos afluentes do rio SNA- 
KE no IDAHO, que é um estado do noroeste ameri- 
cano, numa região onde a água não abunda, já que 
fica na zona semi-desértica entre as montanhas Ro- 
chosas e as montanhas Nevadas, Nota-se à esquerda 
o descarregador, em canal, no extremo inferior do 
qual se pode ver a saída de uma pequena descarga 
do fundo. Também aí se vê, na parte inferior da ou- 
tra margem, o edifício central e o da saída da outra 
descarga de fundo. 
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TÚNEL DE DESVIO, DESCARGA DE 


FUNDO E CONDUTA FORÇADA. (1) — ARGILA SILTOSA COM SEIXOS 


(2) — AREIAS, SEIXOS E CALHAUS COMPACTADOS 


(3) — TOUT-VENANT 
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Fig. 3 
A barragem depois de construída e antes do enchimento da albufeira 


4 A ROTURA 


No dia 3 de Junho, com a albufeira quase cheia, 
observaram-se algumas fontes a jusante da barra- 
gem, deitando água límpida, que não causaram ne- 
nhuma preocupação. No dia 4 de Junho apareceu na 
margem direita junto da parte inferior da barragem 
uma primeira ressurgência com água já turva, que, 
no dia seguinte, já havia aumentado apreciavelmente 
d> caudal (fig. 1). Só por volta das 8,30 da manhã 
é que, efectivamente se dá o alarme porque o caudal 
aumenta, em menos de meia hora, de meio metro 
cúbico por segundo para um metro cúbico por se- 
gundo. Na fig. 4 pode-se ver a ressurgência na parte 
inferior da encosta direita, perto do pé da barragem. 
Entretanto, em cima, por volta das 9 horas, aparece 
uma mancha de humidade no paramento junto à mar- 
gem direita e, pouco depois, começou a correr água 
junto à inserção da barragem com o terreno, forman- 
do-se cerca das 10,15, uma pequena cratera na zona 
da mancha (fig. 5). Portanto, nesta altura existiam 
duas fontes, uma em baixo na rocha e a outra em 
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cima, na barragem, a qual já estava a provocar 
erosão do paramento por volta das 10,30 h da manhã. 
A esta hora, ouviu-se um estampido e formou-se, 
a montante, um vortex na albufeira junto da barra- 
gem, ao mesmo tempo que o caudal escoado aumen- 
tava. Em virtude da rapidez com que estavam a au- 
mentar as infiltrações, os engenheiros do Bureau of 
Reclamation alertaram os xerifes da região e a força 
aérea, e todo o mecanismo de protecção às popula- 
ções começou a funcionar. Deve-se a isso, e com cer- 
teza ao facto de ser de dia, não ter sido este um dos 
maiores desastres da história das barragens. Come- 
caram a ser evacuadas milhares de pessoas e por 
volta das 11h da manhã já estavam desabitadas 
várias povoações inteiras. Só alguns indivíduos que 
estavam a pescar junto do rio é que não foram avi- 
sados e foram surpreendidos pela onda de cheia. 
As 179,45 horas foi decidido encher a cratera com 
enrocamento lançado por dois «bulldozers» que come- 
caram a trabalhar às 11 horas. De qualquer modo, 
a cratera vai aumentando e às 11,30 (fig. 6) os ope- 
radores dos «bulldozers» tiveram que retirar porque 
eles iam escorregando para dentro do fosso. Uns 
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Fig. 4 
Ressurgência na parte inferior da encosta (9,00 h) 


Eriig. à 
Início da cratera na parte superior (10,15 h) 
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Fig. 6 
Cratera quando da queda dos «bulldozers» (11,15 h) 


Fig. 7 
Segunda cratera na parte superior (11,20 h) 
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Aspecto 


Fig. 10 


Fig. di 
da rocha na margem direita após a rotura. Note-se a trincheira corta-águas 
e o aspecto diaclasado da rocha 


ló 


minutos depois forma-se um segundo buraco no pa- 
ramento da barragem (fig. 7), cerca das 11,55 cai 
n coroamento (fig. 8), produzindo-se então o colapso 
total (figs. 9 e 10) que arrastou aproximadamente 
409% do aterro, as obras de jusante (central e câmara 
da controlo da descarga de fundo) deixando à vista, 
por uma larga extensão do rio, a rocha das encostas, 
onde, tal como no local da barragem (fig. 11), se 
pode apreciar o carácter diaclasado e aberto da massa 
rochosa, com predominância das fendas horizontais. 
O ímpeto das águas a grande velocidade, pelo rio 
abaixo, como aliás é corrente nas roturas de barra- 
gens, arrastou tudo, fazendo desaparecer todas as 
obras de juzante (central e câmara de controlo da 
descarga de fundo da M. E.), nada ficando do que 
foi posto pelo homem (foto do frontespício). 


5. EFEITOS DA ROTURA 


Foram inundados 40 000 ha de terras, 780 casas 
foram danificadas e mais de 3000 habitações sofre- 
ram prejuízos, sendo afectadas 4 000 famílias. Se 
não houve muitas mortes é porque, como dissemos, 
toda esta gente já estava avisada havia umas horas. 
Quando a cheia chegou às planícies, havia pelo menos 
umas três horas que andavam por lá helicópteros e 
aviões e que os telefones, polícia e bombeiros, etc. 
actuavam, tudo isto numa movimentação de facto 
impressionante. Aconteceu com esta barragem aquilo 
que aconteceu até hoje em metade dos acidentes 
ocorridos: a albufeira enchia pela primeira vez. 

O acidente saldou-se em catorze mortos, oitenta 
feridos e cerca de trinta milhões de contos de pre- 
juízos. A cheia inundou numerosas estradas e insta- 
lações e afectou cinco povoações. 


Logo que o acidente se deu, foram nomeadas 


(pelas várias entidades oficiais) três comissões de 
inquérito: uma de especialistas de barragens de terra 
e de mecânica dos solos, outra de elementos de orga- 
niemos federais e outra do próprio Bureau of Recla- 
mation. Essas comissões de inquérito consideraram 
já todos os dados e evidências e realizaram cálculos 
e outros estudos, além de ensaios no local, perfura- 
ções, trincheiras de reconhecimento, etc. A primeira 
comissão de inquérito terminou o seu relatório no 
final de Dezembro de 1976, o qual foi imediatamente 
submetido ao Secretário do Interior do Governo Fe- 
deral (Ministro de Fomento) e ao Governador do 
Estado de Idaho. 


6. CAUSAS 


A propósito deste e de outros graves acidentes 
em barragens de terra, alguém já perguntou: «será 
que a rotura de Teton representa uma fraqueza ine- 
rente às barragens de terra”» Em nosso entender 
não há razão nenhuma para responder afirmativa- 
mente. No entanto, estamos certos que, por não 
serem evidentes as causas da rotura, muito se dis- 
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cutirá sobre ela, tanto mais que o Bureau of Recla- 
mation tem dos melhores especialistas mundiais em 
barragens. A rotura não foi devida a nenhum tremor 
do terra. De facto, todos os sismógraios postos na 
barragem e os das quatro estações sismográficos 
da zona, nenhum deles registou qualquer movimento, 
a não ser a rotura da própria barragem. Essa liber- 
tação de energia, sim, está bem registada nos sismó- 
grafos de vários tipos que lá estão, tanto na obra 
como nas redondezas. As investigações já efectuadas 
e outras em curso no sítio talvez possam determinar 
qual foi a linha preferencial de percolação das águas, 
quando essa percolação das águas começou a arras- 
tar material, isto é, a causar erosão interna — o 
chamado «piping>. De qualquer modo, a comissão de 
especialistas entende que as diaclases da rocha e as 
fissuras pre-existentes no aterro, ou a fracturação 
hidráulica deste por pressões da água da albufeira 
vindas das diaclases da rocha, podem ter sido uma 
das causas que originou todo o processo de rotura. 
Essa erosão interna pode ter-se iniciado na inter- 
face entre o aterro e a fundação por erosão sucessiva 
do material impermeável do núcleo, que, dado o seu 
carácter siltoso, é muito erodível, O que não há dú- 
vida nenhuma, pois as fotografias anteriores o de- 
monstram, é que a filtração e a erosão interna co- 
meçam precisamente na zona da trincheira do cor- 
ta-águas da margem direita. Cálculos por elementos 
finitos dos estados de tensão devidos ao peso próprio 
mostram, sem lugar para discussão que, pelo facto 
de existir uma trincheira apertada e profunda em 
rocha rigida, se dá um efeito arco que reduz substan- 
cialmente as tensões do peso próprio da barragem 
no aterro que preenche a trincheira. Tal facto é bem 
conhecido e tem sido causa de problemas graves em 
barragens de terra. A fig. 12, baseada nos cálculos 
da comissão de inquérito, mostra que as tensões no 
aterro dentro duma trincheira apertada são baixas 
e, portanto, nesse aterro, as pressões da água da 
albufeira (que são elevadas) podem provocar fractu- 
ração hidráulica, Aliãs as próprias concentrações de 
tensões nos ângulos podem produzir fissuras. Por ou- 
tro lado, a água circulando nas fendas da rocha em 
contacto com o aterro, como este é siltoso, deve ter 
podido erodi-lo facilmente. Também a dificuldade em 
impermeabilizar, com a injecção, a zona da rocha 
do fundo da trincheira junto da superficie, pode ter 
deixado caminhos, montante-jusante, para a água. 
De qualquer modo, houve várias possibilidades de que 
a água começasse a circular através da fundação 
e do corta-águas com gradientes hidráulicos muito 
fortes que poderiam erodir e destruir o aterro e até 
a própria rocha. 


A experiência mais recente tem mostrado que 
o3 fortes gradientes hidráulicos nos corta-águas ou 
nas fundações das barragens, qualquer que seja o 
seu tipo, a circulação da água na rocha em contacto 
com aterros, a erodibilidade dos núcleos das barra- 
gens de terra, a fissuração desses núcleos (que pode 
ser devida a várias causas, tais como concentração 
de tensões, assentamentos diferenciais, fracturação 
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hidráulica, etc.), a falta de tensões do peso próprio 
(por alívio lateral a que se tem chamado «efeito 
arco» ou por etapas de construção não apropriadas) 
são causas de acidentes muito graves de barragens 


de terra e até da sua rotura. No caso destas obras, 
como tragicamente o provaram em 1976 as barragens 
do Bolan no Paquistão e outras, a insuficiência dos 
descarregadores também é fatal. 


lona mais comprimida 
(efeito arco) 


Zona descomprimida 
(Sujeíita a fracturação 
hidraulica) 


Esquema das trajectória 
das tensoes devidas ao 
peso proprio 


Linhas de igual relação 
entre a tensão vertical 
e a pressao do peso 
proprio por colunas 


(Relatorio da Comissão 
de Inquerito) 


Fig. 12 - ESQUEMA DAS TENSÕES DO PESO PRÓPRIO NUMA DARRAGEM DE TERRA 
FUNDADA EM ROCHA DURA COM UMA TRINCHEIRA DE CORTA-ÁGUAS. 


NOTA DA REDACÇÃO: a redacção da «Técnicas deseja agradecer ao 8r. Professor Laginha Serafim pela prontidão com 
que acedeu ao nosso convite em publicar o texto da palestra «Desastres em barragens em 1976 — o caso Tetons. 
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Feiras, Congressos, Exposições... 


CONGRESSO 77 


A Ordem dos Engenheiros irá realizar, na segunda quinzena de Novembro em Lisboa, um Congresso 
que se destina a tratar de assuntos de carácter técnico, científico e profissional e ainda apreciar a actividade 
desenvolvida pela Ordem. 


Pretendem os organizadores, através deste Congresso, debater os grandes problemas da vida nacional 
ligados à esfera de acção dos engenheiros, contribuindo assim para a resolução desses problemas, visando ao 
mesmo tempo a elevação do nível técnico da engenharia portuguesa e a valorização da sua imagem no 
conceito interno e internacional, 


O Congresso será dividido em duas partes. A primeira, aberta a todos os membros efectivos da Ordem, 
destina-se a apreciar e votar o Relatório e Contas do Conselho Directivo Nacional. A segunda parte, aberta 
a pessoas interessadas e aos membros da Ordem, irá tratar de assuntos de carácter técnico, científico e pro- 
fissional, Nesta parte incluir-se-ão simpósios, conferências, visitas técnicas e ainda outras manifestações de 
carácter técnico e social. 

Os temas do Congresso 77, escolhidos após inquérito à classe, são os seguintes: 


Ti-—-A formação do engenheiro; 

T 2 — Exigências funcionais e processos construtivos para a habitação social; 
T 3 — Unidades de grande potência para a produção de energia eléctrica; 

T 4 — Concepção e projecto de sistemas de telecomunicações no País; 

T 5 — Combustíveis: a escassez mundial e as fontes nacionais; 

T 6 — Metalurgia de metais ferrosos; 

T 7-— Tecnologia de novas culturas; 


T 8 — Aplicação da teoria dos sistemas em engenharia. 


Os temas T1 e 'T 8 serão debatidos em sessões plenárias e os restantes em sessões decorrendo simul- 
taneamente, 

As pessoas interessadas em participar na parte cultural do Congresso (Membros da Ordem, Estudan- 
tes, Outros participantes e Acompanhantes) deverão contactar com: 


ORDEM DOS ENGENHEIROS — Comissão do Congresso 77 
Av, António Augusto de Aguiar, nº 3-D-— LISBOA -1 


4 X BAU FRANKFURT 


De 12 a 20 de Novembro de 1977 realiza-se em Frankfurt, Rep. Federal Alemã, esta feira que incide 
sobre: materiais de construção, elementos de construção, sistemas de construção e acabamentos. 

Para mais informações dirija-se: 

AUSSTELLUNGS — GESSELLSCHAFT MAINZ GmbH & Co. 

K6 — Girtnergasse 16, Postfach 1129, 6500 MAINZ 1, 

REP. FEDERAL ALEMA, 
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S. A. des Ateliers de Sécheron 


1211 Genêéve 21 
Representada em Portugal por 


SOCIEDADE DE ELECTRICIDADE BROWN BOVERI, LDA. 
Rua de Sá da Bandeira, 481-2.º —- PORTO 
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Nas máquinas 

+GF+, de deca- 
pagem ajacto 

por turbina, 


Dispomos da 
melhor técnica 
de decapagem 
para cada cam- 


são deca- Ro Cd A TO ' podeaplicação. 
padosfio máquinaevarãodedimensõese | Existem, também, máquinas +GF+ para a 


pos, capazes de serem trefilados, por pro- chapas, arco, tubos e perfilados incluindo 
cesso económico que defende o ambiente. estruturas metálicas e muito mais ainda. 
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Simulação numérica das estatísticas de Maxwell-Boltzmann, 
Bose-Einstein e Fermi-Dirac, baseada num gerador 
de números pseudoaleatórios DL") 


RESUMO 


Faz-se um estudo comparativo da obtenção por 
simulação numérica, no domínio das pequenas popu- 
lações de particulas, das funções de distribuição cor- 
respondentes às estatísticas de Maxwell-Boltzmann, 
Bose-Einstein e Fermi-Dirac. E uma aplicação dos 
métodos de Monte-Carlo à simulação de um processo 
não determinístico, o das colisões aleatórias de uma 
população de partículas, utilizando um gerador de 
números pseudoaleatórios. Estuda-se a obtenção do 
equilíbrio termodinâmico para populações correspon- 
dentes a uma só estatística, do tipo da lei de radia- 
ção de Planck, e também a interacção entre parti- 
culas de estatísticas diferentes, como electrões e fo- 
nões, 

Este trabalho é um prolongamento de [1]. 


1. INTRODUÇÃO 


A função de distribuição para as três estatísticas 
pode ser sintetizada por: 


n(W) = N(W) [Lexpla + BW,) + AJ-= (1,1) 


onde 
W — energia 


n(W) —- densidade de partículas em energia 

N(W) — densidade de estados quânticos em ener- 
gia 

«5 — constantes a determinar pela energia total 
W, e pelo número de elementos N da população 
d: partículas. 


H. J. S. ABREU SANTOS 
Eng.º Electrotécnico IST 
Assistente do IST 


ABSTRACT 


A comparative study of the obtention of the 
distribution functions corresponding to the Maxwell- 
-Boltzmann, Bose-Einstein and Fermi-Dirac statistics 
is done by numerical simulation on the domain of 
the small populations of particles. It is an applica- 
tion of the Monte-Carlo methods to the simulation of 
a non-determinístic process, that of the random, coli- 
sions of a population of particles, using a pseudo- 
“random number generator. Results are presented 
of the obtention of the eqguilibrium situation for 
populations of particles of the same type, such as 
Planck's radiation law, as well as for different par- 
ticles, like electrons and phonons. 


A — constante que para cada uma das estatiís- 
ticas toma o valor: 


Maxwell-Boltzmann A=Q9 
Bose-Einstein A=—l 
Fermi-Dirac A=1 


Esta diferença, aparentemente trivial, no: resul- 
tado é a consequência das partículas obedecerem 
ou não: 


1) A indiscernibilidade entre partículas 
2) Ao princípio de exclusão de Pauli 


O comportamento das partículas está resumido 
no Quadro I. 


() Publicação do Laboratório de Medidas Eléctricas do IST integrada na Linha 1 do CEUL subsidiada pelo Ins- 


tituto de Alta Cultura. 
(*) Manuscrito recebido para publicação em 14/1/77, 
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QUADRO I 


Estatística Prince, Excl. Pauli 


1) Indiscernibilidade | 2) 
M EB Não Não 
BE |. Sim Não 
“FD Sim Sim 
2. O MODELO 
Leitura Dsdes | 
| Sata DB é 
Entrada AS 
ca 


Ds 
E missao 
Não Y 
Entrada B 7 
Valores med8, 


Não CNENTES 
Sim 
Escrita Result/ 


Fig. 1 


Diagrama de fluxo do modelo de colisão 

N — Contador de colisões 

NMAX — Número de colisões pretendido 

1 — Leitura de dados referentes ao preenchimento inicial 
dos niveis de energia 

&— Escolha do nivel de energia de onde sai a partícula A 

3 — Escolha do nível de energia de onde sai a partícula E 

4 — Determinação da energia da colisão, soma das energias 
dos níveis 

d— Escolha do nível de energia de entrada mais baixa, ao 
qual se atribui a partícula À 

b — Determinação da energia com que fica a partícula E, 
diferença entre a energia da colisão e a energia da 
partícula A após a colisão. Se essa energia for maior 
do que a do nivel mais alto em estudo, a partícula B 
foi emitida e deixou de pertencer à população em 


estudo. 

7— Sé não houve emissão, atribuição da partícula B ao 
nível de cnergia de entrada após a colisão, 

& — Acumulação de valores para determinação do preenchi- 


mento médio e do desvio quadrático médio dos preenchi- 
mentos para todos os niveis. 


mm 20 


Para um número finito de níveis de energia, dis- 
cretos e equidistantes pretende-se obter a distribuição 
de equilíbrio de uma população de partículas, par- 
tindo de uma distribuição inicial qualquer, definida 
pelo número de partículas n, existentes nos N, 
estados quânticos do nível i de energia W.. 


O mecanismo que conduz ao equilibrio é o que 
corresponde ao diagrama de fluxo da Fig. 1, que 
representa uma «colisão» pseudoaleatória entre dois 
elementos da população, com conservação da energia. 
Essa colisão pode ser desdobrada em três partes: 

1) Atribuição pseudoaleatória dos níveis de ener- 
gia donde saiem as duas partículas que vão colidir. 

2) Colisão. 

3) Atribuição pscudoaleatória das partículas aos 
níveis de energia onde entram, com conservação da 
energia da colisão, mas eventualmente distribuida 
doutra forma pelas particulas. 

As diferenças mencionadas em 1. entre o com- 
portamento das partículas pertencentes a estatísticas 
diferentes, traduzem-se no nosso modelo por diferen- 
cas na atribuição dos níveis de saída e de entrada. 


Todas as configurações microscópicas discerni- 
veis que o sistema pode assumir com o número de 
partículas N e a energia W, de que dispõe são equi- 
prováveis. Sempre que se modifica o sistema por 
entrada ou saída de partículas dos ou nos níveis de 
energia, esta modificação tem que respeitar a equi- 
probabilidade das configurações discerníveis. Os es- 
tados quânticos num mesmo nível de energia são 
também equiprováveis. 


214 ESTATÍSTICA DE MAXWELL-BOLTZMANN 
Os blocos referidos em 2.1 são os da Fig. 1 


Bloco 1 


Terão que estar definidos no início para cada um 
dos níveis i (i=1..m) de energia W, os N en. 


Blocos 2 e 3 


sendo as partículas discerníveis, ao retirar uma 
partícula ao sistema posso obter tantas configurações 
diferentes quantas as partículas do sistema. A esco- 
lha da partícula que sai é feita de um modo pseu- 
doaleatório de entre todas as partículas, Assim cada 
nível de energia tem na saida de uma particula um 
peso que é igual ao número de particulas existentes 
nesse nivel. 


Ficam definidos neste bloco as energias dos ni- 
veis de saída Wo € Wan 


Bloco 4 
E calculada a energia da colisão W.o, = Wa + 


+ Won qUe deverá ser a mesma após a colisão, mas 
eventualmente distribuída doutra maneira. 


TECNICA 443 


Blocos 5 e 7 


Uma vez que as partículas não estão submetidas 
ao princípio de exclusão de Pauli, poderão ocupar 
qualquer estado quântico quer já esteja ocupado quer 
não; todas as combinações de pares de níveis que 
respeitem a energia da colisão são possíveis, A esco- 
lha dos níveis de entrada é feita de forma pseudoalea- 
tória de entre todas as combinações de estados quán- 
ticos possíveis. Sendo os estados quânticos dum mes- 
mo nível de energia equiprováveis, o peso de cada 
combinação de niveis de energia é o do produto dos 
estados quânticos do par de níveis que lhe respeita. 


Bloco 6 


Sendo possivel que a energia de colisão seja o 
dobro da energia do nível mais alto em estudo, pode 
a combinação de níveis de entrada escolhida pseu- 
doaleatoriamente corresponder a uma das partículas 
entrar num nível superior ao nível mais alto em es- 
tudo. Haverá emissão de uma partícula com o corres- 
pondente decréscimo de N eW,. 


Bloco B 


Acumulação de dados para cálculo do valor mé- 
dio dos preenchimentos e do desvio quadrático médio, 

Ds blocos 2, 3, 5 e 7 envolvem escolhas pseudoa- 
leatórias de um elemento de uma população dada, 
Para isso recorrem a um gerador de números pseu- 
doaleatórios [2] cuja utilização é fundamental em 
toda a simulação do processo de «colisão». No nosso 
caso, é de lei de repartição uniforme. 


2.2 ESTATÍSTICA DE BOSE-EINSTEIN 
Os blocos referidos em 2.2 são os da Fig. 1 
Bloco 1 


Terão que estar definidos para cada um dos 
níveis i (i=1...m) de energia Wo os n,, N, e as 
ocupações dos próprios estados quânticos dentro de 
cada nível de energia. 

Ao contrário do que acontece com as outras es- 
tatísticas somos obrigados a manter o registo não só 
do número de partículas em cada nível de energia 
mas também do número de partículas em cada estado 
quântico. 


Blocos 2 e 3 


Sendo as partículas indiscerníveis e podendo ha- 
ver mais de uma partícula em cada estado quântico, 
ao retirar uma partícula obteremos configurações 
diferentes desde que retiremos a partícula de estados 
quânticos diferentes. A escolha da partícula que sai 
é feita de modo pseudoaleatório de entre todos os 
estados quânticos que têm uma partícula ou mais. 
O peso de cada um dos níveis de energia é portanto 
proporcional ao número de estados quânticos nesse 
nível com uma ou mais partículas. 
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Blocos 5 e 7 


As partículas não estão submetidas ao princípio 
de exclusão de Pauli, pelo que podem ocupar qual- 
quer estado quântico, quer já esteja ocupado quer 
não, Todas as combinações de pares de níveis que 
respeitam a energia da colisão são possíveis, sendo 
o peso de cada combinação o do produto dos estados 
quânticos do par de níveis que lhe diz respeito, tal 
como nos blocos 5 e 7 de 2.1. 

Temos no entanto que fazer a atribuição das 
particulas na entrada não a níveis de energia mas a 
estados quânticos dentro desses níveis, tendo que se 
manter o registo, utilizado nos blocos 2 e 3, do preen- 
chimento desses estados quânticos. 


Blocos 4, 6 e B 


Idênticos aos correspondentes em 2.1. 
2.3 ESTATÍSTICA DE FERMI-DIRAC 


Os blocos referidos em 2.3 são os da Fig. 1 


Bloco 1 


Terão que estar definidos para todos os níveis 1 
(i=1...m) de energia W, os n, N, e os estados 
quânticos vagos N, —n,. 


i* 


Blocos 2 e 3 


Sendo as partículas indiscerníveis mas não po- 
dendo estar mais de uma em cada estado quântico, 
ao retirar uma partícula obtemos configurações dife- 
rentes retirando particulas de estados quânticos dife- 
rentes, A saída é escolhida pseudoaleatoriamente de 
entre todos os estados quânticos ocupados, necessa- 
riamente com uma única partícula cada. O peso de 
cada nível de energia será portanto proporcional ao 
número de particulas existentes nesse nível. Embora 
por razões diferentes chegamos no entanto ao mes- 
mo resultado que na estatística de Maxwell-Boltz- 
mann, 


Blocos 5 e 7 


Todas as combinações de níveis cujas energias 
somem Wo, Satisfazem. Devido ao princípio de ex- 
clusão de Pauli, as entradas só se podem efectuar 
em estados quânticos vagos. A escolha das entradas 
é feita de forma pseudoaleatória de entre todas as 
combinações de estados quânticos vagos possíveis. 
O peso de cada combinação de pares de níveis de 
energia é proporcional ao produto dos estados quân- 
ticos vagos do par de níveis em causa. 


Blocos 4, 6 8 
São idênticos aos correspondentes em 2.1. 


Os procedimentos anteriores estão resumidos no 
Quadro TI. 

No Apêndice I está demonstrado que a maneira 
d> proceder descrita em 2.1, 2.2 e 2.3 tem como dis- 
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tribuições de equilíbrio, respectivamente, as distribui- 
cões de Maxwell-Boltzmann, Bose-Einstein e Fermi- 
«Dirac. 


QUADRO II 


Saida | Entrada 
proporcional a 


Estatistica 
| proporcional a 


N." de estados 


M B N." de partículas quânticos 


N.” de estados|N.” de estados 


o» di quânticos ocupados quânticos 
N.º de estados 
F D N.º de partículas 


quânticos vagos 


3. RESULTADOS 
3.1 GASES ISOLADOS 
3.1.1 ESTATÍSTICA DE MAXWELL-BOLTZMANN 


Consideraram-se 20 níveis de energia e uma den- 
sidade de estados N,=20, igual para todos os níveis. 
A população era de N =215 partículas que se distri- 


20; 


15 


10 


— o — E o O TES E = — c— — — 


10 0 n; 
Fig, 2 
Maxwell-Holtemann; distribuição inicial 


() O computador utilizado para todos os resultados apresentados foi o do Centro de Cáleulo das Universidades 


de Lisboa, [IBM/360. — 44, 


buiam 100 no nível 1, 100 no nível 4, 14 no nível 5 
e 1 no nível 6, correspondendo-lhes a energia total 
W,=576, tomando para unidade de energia a energia 
que separa dois níveis consecutivos, Ao fim de 20 000 
colisões não tinha havido emissões. 

Na Fig. 2 apresenta-se a situação inicial e na 
Fig. 3 e no Quadro III apresentam-se os resultados 
da média das 10 000 últimas colisões de 20 090 com- 
parados com os valores previstos para N =215 e 
W,=576, com um erro menor que .005, por 1.1, A=0, 
válida no domínio dos grandes números, tendo q e f 
sido determinados por 3.1 e 3.2. 

N dy N, cexpía =8W,) a (3.1) 


u 


i 


Os tempos de CPU do programa foram aproxi- 
madamente de 1 minuto para 20000 colisões (1). 


W 


k 


20. 


15 


10 


1 
| à á 
Fig. 3 
Maxwell-Boltzmann; distribuição final 
Prevista 


O 20000 colisões; média de 10 000 
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QUADRO III 


Nivel Previsto 10 000 col 
1 80.21 TT.8l 
2 50.29 51.46 
3 31.53 32.29 
4 19.77 20.76 
3) 12.39 12.05 
6 T.iT 8.58 
T 4.87 4.59 
8 3.05 3.02 
9 1.02 1.85 

10 1.20 7 
11 15 33 
12 47 al 
13 30 so 
14 19 AB 
15 AZ 18 
16 OT 09 
11 05 06 
18 03 DO 
19 02 09 
20 01 DO 


3.1.2 ESTATÍSTICA DE BOSE-EINSTEIN 


Consideraram-se 100 níveis de energia, e uma 
função densidade de estados em energia, N, corres- 
pondente à parábola N =12(W, —.5)º de eixo de 
simetria nos N.,. 


erro apreciável, São muito pouco prováveis as situa- 
ções com mais de uma partícula num mesmo estado 
quântico. 


Dum ponto de vista estatístico é irrelevante o 
que são os n, N, e W.. 


Se n, for o número de osciladores (corpuscular- 
mente, fotões) para a frequência própria f, corres- 
pondente à energia W,=hf,, a distribuição em equi- 
líbrio dos n, pelos W, conduz, numa representação 
de n/W, em função dos À, =hc/W, a uma distribuição 
da densidade de energia do tipo da Lei de Radiação 
ds Planck. 


Difere desta por termos forçado um número de 
partículas N constante, o que não é verdade para 
fotões. Assim, (1.1) para fotões teria ,=0. 


No modelo utilizado atingiu-se o equilíbrio por 
trocas de energia entre partículas. Para fotões, dada 
a linearidade das equações de Maxwell, estes não 
interactuam, só sendo possível atingir-se o equilíbrio 
na presença de matéria [3]. 


A população foi de N = 156 partículas, não tendo 
havido emissões, e a energia W =7 51. A situação 
inicial, em termos dos preenchimentos dos diversos 
níveis de energia, e o número de estados quânticos 
por nível de energia foram as do Quadro IV. Acima 
do nível 8 não havia partículas, e tirando o nível 1, 
as partículas ocupavam todas estados quânticos dife- 
rentes. O comportamento inicial será portanto um 
comportamento à Maxwell-Boltzmann. 


QUADRO IV 


W,=751 W,=915 
Nível Estados Electrões Estados ocupados Electrões Estados ocupados 
1 Ni o! J 4 ho) 
2 27 1 1 17 17 
3 T5 1 1 21 21 
4 147 20 20 aa 22 
5 243 25 25 20 20 
6 363 24 24 17 17 
7 507 24 24 14 14 
8 675 24 24 11 11 
9 867 O) LO) 8 8 
10 1083 o Li) 6 6 
al 1323 0 O 3 3 
12 1587 0 O O) O 
13 1875 0 O 8 B 
14 2187 0 É) 4 4 


Para que o programa fosse praticável em tempo 
de CPT, os níveis 11 a 100 foram considerados como 
se se tratasse de níveis de Maxwell-Boltzmann. Como 
as duas estatísticas coincidem quando em cada estado 
quântico existe só uma partícula e no exemplo esco- 
lhido, no nível mais desfavorável, nível 11, são pre- 
vistas em equilíbrio menos de 5 partículas em 1323 
estados quânticos equiprováveis, não haverá por isso 
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Os valores médios das 19000 últimas colisões 
de 25 000 estão comparados no Quadro V e na Fig. 4 
com os valores previstos por 1.1, A=-1, com q e B 
calculados por 3.3 e 3.4 com um erro menor que .005. 


| 
No = D Ni [exp(r+ 8W,) 17? 


(3.3) 


| 
Wo= ON W, [exp(a+8W))—1] (3.4) 


23 


100% 


20 


& 
a 
o Fig. 4 


z Bose-Einstein; distribuição final 
(2) — Prevista, W,=751 x 
(1) 


Prevista; W,=915 O 
E ea de pEM Prevista; W,=751; Maxwell-Boltzmann 


()— — —— — — Prevista W,=-7bl; N, —l 


25 000 colisões (W,=7 51) 
20 000 colisões (W,=915) 


D! mad 


QUADRO V 


Nivel Previsto 25 000 col, Previsto 20 000 col. 
(med. 10 000) (méd. 15 000) 
1 13.88 11.89 4.85 4.99 
2 2217 22.12 15.67 16.12 
3 24.63 26.94 20.95 20.83 
4 23.05 23.16 21.93 22.10 
5) 19.51 18.86 20.32 19,84 
6 15.42 16.33 17.45 18.12 
fi 11.59 11.75 14.22 14.55 
8 8.37 7.84 11.13 10.99 
9 5.86 6.00 8.44 8.21 
10 4.00 3.82 6.24 6.68 
11 2.68 247 4.02 4.65 
12 1.76 1.64 3.22 3.11 
13 1.14 1.26 2.26 2.13 
14 3 8 1.57 1,44 
15 46 42 1.07 1.02 
16 .29 33 3 60 
LT 18 AO 49 57 
18 11 Rio) So 25 
19 OT OS 22 26 
20 D4 05 14 15 
21 OS DO 09 12 
22 OZ 02 06 04 
23 O “0 Od O 
24 01 Ri) 03 01 
25 OO 00 02 O 
26 00 09 0 01 
27 DO 00 01 01 
28 00 00 DO 00 
29 00 OO DO 00 
So DO DO OO 00 


Encontra-se na mesma figura a distribuição de 
equilíbrio prevista, N —156 e W,='751, se o compor- 
tamento das partículas fosse de Maxwell-Boltzmann. 

Numa segunda situação, com o mesmo número 
de partículas mas W =915, a situação inicial foi a 
do Quadro IV, com as partículas em estados quân- 
ticos diferentes excepto no nível 1, e sem partículas 
acima do nível 14, 

Os valores médios das 15000 últimas colisões 
de 20000 encontram-se no Quadro V e na Fig. 4 
comparados com os previstos por 1.1 (A=-—1l) 3.3 
e 3.4 com um erro menor que .005. 

Os tempos de CPU foram aproximadamente de 
& minutos para 20 090 colisões. 

Os valores obtidos para pequenas populações es- 
tão próximos dos previstos por expressões aplicáveis 
para grandes populações. Não é no entanto seguro 
que estas expressões sejam válidas no domínio dos 
pequenos números. Em relação ao mecanismo usado 
nas colisões, a demonstração feita no Apêndice I de 
que tem como distribuição de equilíbrio a de Bose- 
“Einstein, pressupõe também que N +n, —1 = N,+ 
+n, para determinação da probabilidade de saída. 
E portanto possível que o mecanismo de simulação 
usado conduza para pequenas populações a uma dis- 
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tribuição diferente, Na impossibilidade de trabalhar 
com n, maior nos baixos níveis, porque dado o seu 
crescimento parabólico tomaria nos níveis mais altos 
valores que obrigariam no programa a uma ocupa- 
ção de memória demasiado grande, fizeram-se ensaios 
a N, constante igual para todos os níveis mas com 
um valor superior ao usado no exemplo anterior para 
os níveis mais baixos, tendo-se verificado um melhor 
acordo entre valores encontrados e previstos para 
grandes populações. Por outro lado tentou-se uma 
correcção dos valores previstos trabalhando por simu- 
lação com N,, mas fazendo a previsão para a mesma 
população e energia mas como se cada nível tivesse 
N,—l estados. (Fig. & curva (2) ). Esta correcção 
nunca poderia ser inteiramente correcta porque as 
saídas são realmente previstas para N,—l, mas as 
entradas são também previstas para N—l e não 
para N.. A correcção mostra-se no sentido correcto, 
embora insuficiente, 


Repare-se que todas as colisões em que as saídas 
correspondam aos pares de níveis 1;1 ou 2;1 têm 
como níveis de entrada os de saída para que haja 
conservação da energia, permanecendo o sistema com 
preenchimentos dos níveis de energia iguais aos 
anteriores à colisão. Num caso extremo em que só 
houvesse partículas nos níveis 1 e 2 o sistema per- 
maneceria sempre com a distribuição inicial, que a 
não ser a de equilíbrio, nunca permitiria que este se 
verificasse. Este aspecto comum a todas as estatisti- 
cas, é tanto mais relevante quanto, dado o pequeno 
número de níveis, as probabilidades de intervenção 
dos níveis em colisão variam muito de nível para 
nível. Se forem maiores nos níveis mais baixos, 
grande parte das colisões mantêm o sistema inva- 
riante. 


3.1.3 ESTATÍSTICA DE FERMI-DIRAC [1] 


Consideraram-se 20 níveis de energia e 20 esta- 
dos quânticos por nível. A situação inicial (Fig. 5) 
correspondia a todos os estados quânticos ocupados 
até ao nível 9 inclusivé, 10 estados ocupados nos 
níveis 10 a 16, e 5 estados no nível 17 (N,=255; 
W,=1895). Na Fig. 6 e Quadro VI apresentam-se 
os valores média das 10000 colisões entre as 
10000 e 50000, em que a população se mostrava 
estável após sucessivas emissões, 243 electrões com 
uma energia de 1625, comparados com os resultados 
previstos por 1.1, A=1, com q e 8 calculados por 3.5 
e 3.6 com um erro menor que .005. O W e oN, 
utilizados foram os da situação final, 


i 
No= > N, [exp(a+ BW, )+1)" (3.5) 


N 


<a 


i 
= DD NWi[ep(«+8Wi)+1] (39) 


25 


Fig. 5 
Fermi-Dirac; distribuição inicial 
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ta 


a” 


] 
— ij —— dj == A a 


o q 5 10 15 20 Nimed 
Fig. 6 
Fermi-Dirac; distribuição final 
Prevista 
e 50 000 colisões 
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QUADRO VI 


Niver Previsto 50 000 col. 
1 20.00 20.00 
- 20.00 20.00 
E 20.09 20.00 
4 20.00 20.00 
5 20.00 20.00 
6 19.99 19.99 
7 19.96 19.98 
8 19.89 19.87 
9 19.65 19.64 

10 19,90 18.98 
11 17.24 17.19 
12 13.46 15.22 
13 8.07 7.89 
14 3.64 4.3s 
ló 1.36 1.05 
16 47 49 
17 16 17 
18 O OO 
19 02 OO 
20 01 OS 


3.2 INTERACÇÃO DE PARTÍCULAS OBEDECENDO 
A ESTATÍSTICAS DIFERENTES: «BANHO» DE 
BOSE-EINSTEIN E GAS DE FERMI-DIRAC 


O modelo utilizado foi descrito em pormenor 
em [4]. Basicamente corresponde a colisões entre 
pares de partículas, uma do gás e outra do «banho» 
à temperatura T, suposto com uma capacidade calo- 
rífica infinita. Completada a colisão, a energia do 
par de partículas após a colisão é a mesma que antes 
mas eventualmente repartida doutra maneira. Feito 
o regresso das partículas, uma ao gás e outra ao 
banho, este repõe a partícula que o abandonou na 
situação energética que tinha anteriormente. Sendo 
assim, a distribuição do banho mantém-se invariante 
com preenchimentos iguais aos valores médios corres- 
pondentes ao equilíbrio à temperatura T. O equilíbrio 
é atingido quando o gás estiver à temperatura do 
banho. 

Seja o gás A de Fermi-Dirac, e o gás E, banho, 
de Bose-Einstein. Cada colisão envolve uma partícula 
do gás A e uma partícula do gás B. Na saída do 
gás A, o peso de cada nível é o do número de parti- 
culas nesse nível. Na saída do gás B, o peso de cada 
nível será o do número de estados quânticos com uma 
ou mais partículas. É portanto necessário conhecer 
à priori, em equilíbrio à temperatura T, qual a pro- 
babilidade de para cada nível encontrar, nos N, 
estados quânticos, estados com uma ou mais partí- 
culas. Pode ser calculada (ver Apêndice II) pela ex- 
pressão n,N,/(n, + N, —l). Basta conhecer os N e 
os n, de equilíbrio à temperatura T pretendida para 
o banho, 

O exemplo apresentado corresponde a uma situa- 
ção bidimensional em que as densidades de estado 
em energia, dada a quantificação no espaço dos mo- 
mentos e as relações entre energia e momento no 
caso de electrões (W=p'/2m com p-— momento 
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linear; m— massa eficaz) e no caso de fonões 
(W=up com u-velocidade do som), são respectiva- 
mente N,, constante e N,, crescendo linearmente 
com a energia. 

Consideraram-se 100 níveis de energia, N,,=29 
e N,,-10 W, sendo W, expresso em intervalos de 
energia entre níveis consecutivos. Os elementos da 
população de Fermi-Dirac eram N =204 e encon- 
travam-se inicialmente na situação da Fig. 7. A po- 
pulação de Bose-Einstein estava em equilíbrio à tem- 
peratura T, tal que a energia característica dessa 
temperatura, KT, era de 5 intervalos de energia. 


Wi 
100 
99 


20 


15 


10 


1 5 O 15º niinic 


Fig. 7 
Gás de Fermi-Dirac em contacto com «banho» de Bose- 
“Einstein, Distribuição inicial 
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resultados das médias das 10000 últimas colisões 
de 50 000 e das 30 009 últimas de 70 000 comparados 
com os valores previstos por 1.1, A=1, com 8=1/kKT com um erro menor que .005. Na Fig. BA apresen- 


calculado a partir da temperatura T do banho e q à ta-se a distribuição em energia para as partículas do 
partir de «banho». | 


Na Fig. 8 e no Quadro VII apresentam-se os N= > N[exp(«+8W)+1]" 7 


Wi 
100: 


ay 


Fig. 8 
| é Gás de Fermi-Dirac em contacto com «banho» de Bose- 
o “Einstein. Distribuição inicial 
2 5 Previsão 
» DO 000 colisões 
x 70 000 colisões 


4 5 0 5 17n 


it med 
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Fig. BA 
Gás «banho» de Bose-Einateln 
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QUADRO VII 


Nível Previsto 50 000 col, TO 000 col. 
(méd. 10 000) (méd. 30 000) 

1 17.16 16.80 16.96 

2 16.64 16.60 16.60 

3 16.04 16.10 16.02 

4 15.37 15.63 15.45 

) 14,62 14.67 14.52 

6 13.80 13.90 13.89 

7 12.91 13.17 12.92 

8 11.87 11.40 11.78 

9 11.00 10.64 10.84 
10 10.00 10.18 10.11 
1 8,09 9.00 9,05 
12 8.03 8.59 8.27 
13 7.09 6.53 6.71 
14 6.20 6.40 6.37 
15 5.38 Del o.26 
16 4.63 4.21 4.68 
17 3.96 3.86 3.99 
18 S.Sb 2.94 2.97 
19 2.84 Aki 2.80 
20 2.38 2.41 2.43 
21 2.00 1.55 1.80 
22 1.66 1.85 1.80 
23 1.38 1.80 1.58 
24 1.15 1.39 1.23 
25 DO To T6 
26 TB .B5 TT 
27 65 96 93 
28 oa Ri 53 
29 44 Rsja. 48 
30) 36 SO 46 
31 80 43 39 
32 24 19 25 
RS] 29 46 So 
Sd 16 Ao 12 
ao 13 09 0s 
36 11 27 20 
31 09 17 14 
3a 07 2 01 
39 06 05 06 
40 OS 11 10 


Acima do nível 4) o preenchimento previsto é 
menor que .05 e nos 60 níveis restantes só dois dos 
preenchimentos médios atingiram valores desta or- 
dem de grandeza, ambos para a média de 10000 
colisões, pelo que não são apresentados valores para 
estes níveis. 

Na Fig. 8 e no Quadro VIII apresenta-se o des- 
vio quadrático médio c* comparado com o valor 
previsto por [5]. 

E] 


E EM as (1— D,.q/N () 


(3.8) 


em que n,., é o número de electrões no nível i em 
equilíbrio termodinâmico. Tal como previsto as flu- 
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tuações são maiores na vizinhança do nível de Fermi 
(nível 10). Acima do nível 40 os valores previstos 
eram menores que .05 e nos valores encontrados nos 
restantes 60 níveis só um atingia este valor, tanto 
para 50000 como para 70000 colisões, pelo que os 
valores também não são apresentados. 


QUADRO VIII 


Nivel Previsto 50 000 col. 70 000 col. 
(méd. 10 000) (méd. 30 000) 

1 2.43 2.62 2.66 

2 2.80 2.26 2.41 

3 3.17 2.36 2.81 

4 3.56 3.07 Sls 

õ 3.93 ST 3.80 

6 4.28 4.38 4.18 

7 4.58 3.66 3.81 

B 4,81 4.18 4,37 

9 4.95 4.28 4.97 
10 5.00 4.76 4.90 
11 4.95 4.91 4.87 
12 4,81 4,85 5.41 
13 4.58 3.85 4.14 
14 4.28 4,05 a. TO 
15 3.93 414 3.98 
16 3.06 2.73 3.12 
17 3.17 2.97 3.20 
18 2.80 2.19 2.19 
19 2.43 2.50 2.44 
20 210 2.83 2.69 
21 1.80 1.04 1.59 
22 1.53 1.61 1.69 
23 1.29 1.30 1.33 
24 1.08 1.23 1.06 
2 SO Nf: 1 
26 o BS To 
27 63 6 1.38 
28 a 4 66 
29 43 2 42 
30 SO 37 38 
31 29 .26 28 
32 24 21 30 
E) 20 4o 37 
Sd 16 09 JO 
JO A3 10 09 
30 Ji 21 Ar 
37 09 16 Lá 
38 97 OZ Rio! 
39 06 03 06 
40 05 AO 09 


O tempo de CPU foi de cerca de 6 minutos para 
TO 000 colisões. 

Os valores encontrados apresentam um bom acor- 
do com os valores previstos para grandes populações, 
melhorando se tomarmos médias sobre um número 
crescente de colisões. 
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Fig. 9 


Gás de Fermi-Dirac em contacto com «banho» de Bose-Einestein. 
Desvio quadrático médio 
Previsto 
o 50 000 colisões 
70 000 colisões 
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4. CONCLUSÕES 


Para qualquer das estatísticas foi bom o acordo 
entre os resultados obtidos para pequenas populações 
e os resultados previstos por expressões deduzidas 
no domínio das grandes populações, embora esse 
acordo não seja perfeito. O acordo é tanto melhor 
quanto maior for o número de colisões sobre o qual 
a média é tomada. O método não é no entanto prati- 
cável para populações da ordem de grandeza das dos 
electrões nos cristais, visto que com os meios de 
cálculo actualmente disponíveis os tempos de CPU 
são da ordem dos ms por colisão e o número de 
colisões para tomar valores médios, se considerarmos 
que é aceitável que cada uma das partículas tenha 
probabilidade diferente de zero de encontrar todas 
as outras, será da ordem do quadrado do número 
de elementos da população. 

Este método de Monte Carlo pode no entanto ser 
útil em situações em que as populações em causa 
sejam pequenas e em que portanto as expressões 1.1 
não são necessariamente válidas. 

O método usado esclarece quais as diferenças 
profundas de comportamento entre partículas de 
estatísticas diferentes, e que não são aparentes nas 
diferenças triviais entre as expressões que traduzem 
os resultados. Em particular, comparando as esta- 
tísticas de Maxwell-Boltzmann e Bose-Einstein a dife- 
rença subtil entre as particulas serem respectiva- 
mente discerníveis e não discerníveis, traduziu-se 
afinal em que na estatística de Maxwell-Boltzmann 
obtêm-se configurações microscópicas distintas reti- 
rando qualquer partícula, mesmo que haja várias 
partículas ocupando o mesmo estado quântico, en- 
quanto na estatística de Bose-Einstein só se obtêm 
configurações distintas retirando partículas de esta- 
dos quânticos diferentes. 
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APÊNDICE | 


Consideraram-se níveis de energia de i=1..m. 
Sejam para cada nível de energia 


W, — energia no nível i 

N, — número de estados quânticos no nível í 
n, — número de partículas no nível i 

N — número total de partículas 


o 


Para o modelo de colisões descrito em 2. uma 
distribuição será de equilíbrio quando, para cada um 
dos níveis, a probabilidade de sair uma partícula desse 
nível for igual à probabilidade de entrar uma par- 
tícula nesse nível após a colisão (1). 


Não vamos demonstrar que os mecanismos des- 
critos em 2.1, 2.2 e 2.3 conduzem às respectivas dis- 
tribuições de equilíbrio, mas só que essas distribuições 
são de equilíbrio para os mecanismos usados, embora 
não sejam necessariamente as únicas. 


A.1 — Maxwell-Boltzmann 


A probabilidade de saída de uma partícula do 
nível i é 


Ps;,-C,n, (41,1) 


na 


onde C, — 1 / > n=1L/N, 


j=l 


Para que haja possibilidade de entrada de uma 
partícula no nível i é preciso que a energia da coli- 
são W.., > W;,- Tendo as partículas em colisão saído 
dos níveis i e k, Wcol=W, + W, que para igual 
afastamento entre níveis de energia é também 
Was — Wi A entrada, após a colisão, ao nível i 


corresponde um peso N,N,, visto que, havendo con- 

servação da energia, se uma partícula entra no nível 

ia outra entra no nível k. São possíveis as entradas 
461 

correspondentes a 3 Nu Nicky. A probabilidade 
n=1 


de a energia de colisão ser W, E é 


Il it-t 


g 
>, Psn Psitkn > Cn, Njtkem 


n=1 == 


Como k pode variar de um até m temos para probabilidade de entrada no nível i 


Po DO ES (41.2) 


Sendo para uma distribuição de Maxwell-Boltzmann 
n,=N, exp (a + AW,) (A1.3) 
de (Al.1) e (A1.2), atendendo a (A1.3) ea que Wi, W,tW, 


psi=CN exp (a + AW,) 


E+-k-l 


m 2 N n Nipken EXP má im | W:)] Ny Ny 
. n==1 E 
egg 2 El 
k==1 
2 N a Dirk nu 


n=1 


1+k-l 
b2 Ny Nigk-n 
— CN exp(a BW, JB mep(e+ tu) Fi — (A14) 
b3 N a Nisk-a 
n=1 


m 
Sendo +. N, exp ( a +BWi ) = Mn = 1/6, 
k=1 


Pp, = (1/N9)N, exp (a + BW) = Ps; 


A.2 — Bose-Einstein 


A probabilidade de saida é proporcional ao número de estados, dos N, estados quânticos, que tiverem 
uma partícula ou mais e é (ver Apêndice II) 


n,N;/(N, + a, —1 ) 


ou, se N + n,d1 


Pq = € 2,N/(N +20) (A2.1) 


onde 


TD» | nn; [ (mA + N; | (A2.2) 


i=1 


A probabilidade de entrada, por um raciocinio análogo ao seguido em A.l é 


I-Hk-1 
RL) r ra 
= > Ca Ma No Pipken Nikon Ni Ny 
P = >» - RR re sat (A2.3) 
k=1, 


(no + N a) (Miro + N Lts) > No à PP 
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Sendo para uma distribuição de Bose-Einstein 


n,=N,/[exp (a + AW,) —1] (AZ4) 


substituindo (A2.4) em (A21) e (A2.2) 
Pa =C N, exp(— «— BW,) 
m 
C=1/ N; exp(-2-BW, ) 
i=1 


Substituindo (A2.4) em (A2.3) e simplificando 


i+ k-d 


o Na Nipin exp | —20— 8(Wi+ Wi) | NM, 
a n=1 
Fpp > ill | 
k=1 É 
a +, Ni k-m 
n=1 
L4-k-1 
N. Ni 
m 2 né“ d+ km 
4 r - h - E. 4 Tl= 
= Cc N, exp re a—BW,) > Ny exp |[— a — BW) a 
k=1 , 
N, Ne k-n 
n=1 


sendo >, Ny eXp (— a— BW, ) = LC; 
k=1 


Pp — C, N, exp(— a — BW,) = Pai 


A.3 — Fermi-Dirac 


A probabilidade de saída é 
Ps, =€Cn, 


onde 
C = 1/N, 


a 


As entradas no nível i só se podem dar nos estados quânticos vagos N, —n,. Tendo isto em conta, por 
um raciocínio semelhante ao seguido em Al, a probabilidade de entrada é 


i+k—l 

q | Pad 4 Ea 4 

ii b> Do Djs km UN = My) (Ny = n,) 
>. tra n=1 
Puto DA Ta 


E (um sata 


n=1 


bd Dj Ms km 


" | hi] Fa bp | 
= 2 tai ST TT a 
k=1 


3 Na 7] ) No, pn | 
nm cr = 
un A--Esm no ; NiLk-n 


n=1 
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sendo para uma distribuição de Fermi-Dirac 
n,/N, = 1/[expla + AW,) + 1] 


então 
n,/N,—1 = exp(a + BW.) (A3.2) 


Substituindo (A3.2) em (A3.1) e simplificando 


itk—l 
um 2 Cu Pitkoa 
Po = Cj (N,-my) ba (N, - He ) E ERR] 
pera exp [2 a B (W, + Wi) | A Da Dicio 
n=1 


Atendendo a (A3.2) 


mm 


ba No (Nf ne—l ) exp (— a —B Wi)=N=1 / CG. 
k=1 


Pr; = CD, = Ps; 


“oe H 
(+) 


Porque uma partícula pode abandonar um nível quer como primeiro quer como segundo termo da coli- 
são, e pode entrar num nível como primeira ou segunda partícula após a colisão, sendo as duas partículas 
indiscerníveis quanto à ordem porque saem ou entram, mesmo no caso de Maxwell-Boltzmann, basta para 
estudar o mecanismo de colisão considerar a entrada e a saída de uma das partículas. Só não ficam contem- 
plados na demonstração os casos em que ag partículas saem ou entram as duas do ou no mesmo nível, Consi- 
derando essas hipóteses teriamos: 


Maxwell-Boltzmann 


Saida 
Ps, = Ps; * Pg; (0,—D/MN —D) =Ps, sen «N, 
Entrada 


Mantém-se 
Bose-Einstein 


Saída 


; es IN 
Ps= Ps + Ps ES rei | a 
: | a RL nm N, 
Nittap=l Jo; é ci 1 
L Ip> No +n, 


se P,, “1, visto que o 2.º termo é da ordem de p,, 


Entrada 
Substituir as hipóteses de entrada 


|rk—l 


> N np Nipt- [o 
n=1 
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por 


i+k—l itk=l 
> Na Nik-n Te Na) 12 [ NH) fg + 1 I= x No Nisy-a se Neã) a» I 
n=1 n=1 
n=5=(1+k)/2 
Fermi-Dirac 


Saída 
Como em Maxwell-Boltzmann 
Entrada 


Substituir as hipóteses de entrada 


I+k—1 
2 (Na == ) (N, +k-n Mpk-n) 
n=1 
por 
i+k— | 
; | 
(Na — Ng) (Nin — Mpa-n) + [ Nua) 2 — May) BM pp — My —1 |] 
n=1 
n&(i-+k)/2 


que é aproximadamente igual ao anterior se as vagas no nivel (i + k)/2 forem muito maiores do que um, 


APÊNDICE II 
Seja um nível genérico de energia com N estados equiprováveis e n partículas indiscerníveis, Cada estado 
pode ser ocupado por uma ou mais partículas. 
Seja k o número de estados com uma partícula ou mais; N-k estados estão vazios. Para garantirmos 
que os k estados têm pelo menos uma partícula, imobilizemos uma em cada um deles, As n-k partículas res- 


tantes podem ser dispostas nos k estados dum número de maneiras que corresponde às combinações com repe- 
tição de k, n-k a n-k 


n=-1 
*n=-k 


Com N estados podemos formar um número de grupos de k estados diferentes 
E N 
k 
O número total de estados com zero partículas será então 
nl AN 
SN) Co li 
k 
Sendo o número total de estados 


cN+a- 


n 


a probabilidade de obter estados vazios é 
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Z(nN)= 


Simplificando e atendendo a que 


B Cha E a 


então Z (1,N) = (N—1)/(N +n—1) 


Ren dE 


e portanto a probabilidade de um dos N estados ter uma ou mais partículas é 


1— Z (0,N) =n0/(N +n—1) 


Sumários de números anteriores 
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1 a 420 (inclusive) — 
421 a 436 (inclusive) — 
438, 439, 442 — 


10$00 cada exemplar 
35800 cada exemplar 
60$00 cada exemplar 


440/441 — 100800 cada exemplar 


TÉCNICA 305 — Outubro 1960 


Valorização individual e progresso colectivo .. 

Bases teóricas para a determinação experimental 
do espectro de potência 

Problemas do estado plástico ... 

Reacções químicas sucessivas de péimeira: ordem 
— constantes cinéticas ... .... 

Cálculo a rotura de estruturas nspacogGnaas pe 
trópicas 


TÉCNICA 306 — Novembro 1960 


Semana de recepção aos novos alunos ... 

Considerações sobre a determinação da carga 
térmica dos motores de combustão interna 

O controlo da fabricação do betão e a medição 
da sua qualidade no estaleiro . 

Amplificação de sinais eléctricos prece espia 
Estudo da amplificação distribuida ... ..... 

Cálculo à rotura de estruturas heterogéneas orto- 
trópicas 


Luis Almeida Alves 


M. J. Abreu Faro 
M. F. Almeida e Sd 


A. G. Silva Pinto 


W. Olezak e A. Sawczuk 


Luís Almeida Alves 

F. de Almeida e Castro 
4, de Sousa Coutinho 
J. M. Jervis Pereira 


W. Olezak e A. Sawczuk 


TÉCNICA 443 


NUMERO 443 OUTUBRO DE 1977 


ANO Lil 


VOLUME XXXIX 


C. D. U. 552.1.067.1 


Ocorrência de scheelite em rochas calco-silicatadas 
da região de Vila Franca das Naves 


RESUMO 


A presente nota refere a ocorrência de scheelite 
em bancadas de rochas calco-silicatadas das zonas 
de Rocamondo e de Freches (carta geológica 18-A). 
Estas rochas constituem um conjunto grosseiramente 
paralelo com direcções N70"W verticais a sub-verti- 
cais. 

Mineralogicamente trata-se de rochas com tex- 
tura granoblástica a porfiroblástica essencialmente 
constituídas por plagioclase, clhinopiroxena diopsídica 
e anfibola, a que se associam, idocrase e granada. 
Acessoriamente estão presentes esfena, goizite, clino- 
goizite, apatite e opacos. A scheelite ocorre em quase 
todas as rochas em pequenos grãos. 

A formação e evolução das rochas calco-silica- 
tadas estará relacionada com acções de metamorfis- 
mo de contacto exercidas pelo granito no seio do 
qual se localizam os retalhos xistosos que englobam 
os níveis estudados. 

A mineralização scheelítica ocorre normalmente 
disseminada no seio da rocha escarnítica, parecendo 
tratar-se de ocorrências pobres do ponto de vista 
econômico não só pela irregularidade da mineraliza- 
ção como também pelo seu grau de finura. 


| — INTRODUÇÃO 


A presente nota inicia um trabalho mais amplo 
e inscreve-se em futuros estudos de pormenor que os 
autores pretendem levar a cabo relativamente a al- 
gumas dezenas de afloramentos de rochas calco-sili- 
catadas mineralizadas por scheelite. 

A scheelite dos afloramentos referidos neste tra- 
balho foi pela primeira vez detectada em trabalho 
de campo efectuado por um dos autores (A, F, Bar- 
ros) e M. M. Azevedo, 
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SUMMARY 


The calc-silicate rocks are concordant with the 
schists which are orientated N 70ºW vertical to sub- 
-vertical, 

These rocks have a granoblastic to porphyro- 
blastic texture and they contain plagioclase, diopsidio 
clinopyroxene, amphibole and some idocrase and gar- 
net. The following accessory minerals were also 
found: sphene, goíisite, clinozoisite, apatite and opacs. 
Scheelite was found in small grains in most of the 
rocks stiwlied. 

Some big schist xenoliths occur inside the gra- 
nite of the area and it is thought that the scheelite 
mineralization and the last part of the evolution of 
the cale-silicate rocks is related to the effects of the 
granite contact metamorphism. 

The scheelite generalhy oceurs in small scattered 
grains in skarns, so it seems that the economic in- 
terest of this mineralization will be small. 


|l— GEOLOGIA E MINERALIZAÇÃO 


A região de Vila Franca das Neves cuja carto- 
grafia geológica está feita à escala 1/50 000 (carta 
geológica 18-A) é fundamentalmente ocupada por 
rochas de natureza granítica e alguns retalhos de 
rochas xistentas que se admite possam pertencer ao 
complexo xisto-grauváquico. (Teixeira, C., Medeiros, 
A, €., Pilar, L., Carvalhosa, A. e Ferro M,, 1963). 

Em duas dessas manchas (Rocamondo e F're- 
ches) ocorre scheelite associada a tipos litológicos 
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calco-silicatados. Estas rochas, mais perfeitamente re- 
presentadas nos afloramentos de Rocamondo, com 
possanças variáveis desde poucos centímetros até 3 
metros ocorrem em camadas intercaladas nos xistos. 
A rocha granitóide de grão grosseiro que envolve os 
retalhos xistentos penetra-os profundamente sob a 
forma de apófises que podem cortar as faixas calco- 
-magnesianas, 

As bancadas calco-silicatadas constituem um con- 
junto grosseiramente paralelo de direcção próxima 
de N70ºW e de pendor vertical ou sub-vertical, con- 
cordante com a rocha intercalante (Fig. 1). 


No entanto e apenas na zona de Rocamondo, 
dada a extensão do afloramento calco-silicatado e a 
presenca, embora em teores baixos, de mineralização 
scheelítica em vários locais, é de pensar ter interesse 
a execução de clássicos trabalhos de pesquisa. 

Esta opinião é reforçada pela existência neste 
local de concessões de volframite, actualmente na 
situação de abandono. 

Do ponto de vista de beneficiação, considera-se 
o minério do tipo silicioso com ocorrência scheelítica 
fina, (A. M. Cerveira, 1971), 


+ + + 


A 4 ar E Granito porfiróide de grão grosseiro 


1! ! Granito de grao muito grosseiro 


T T 7 frequentemente porfiróide 


144! Complexo xisto - granitico - 
5 s 4 -migmatitico 


=” Rochas calco-silicatadas 


0] Concessões mineiras abandonadas 


Fig. 1 - Esboco geológico da região de Rocamondo 


O seu modo de ocorrência é variável, sob a forma 
de finas camadas centimétricas intercaladas no mi- 
caxisto numa extensão de vários metros, ou apresen- 
tando-se em camadas espessas que revelam variações 
composicionais. 

Na mancha de Freches os tipos calco-magnesia- 
nos estão mais pobremente representados e as cama- 
das dificilmente se conseguem observar no terreno, 
sendo no entanto abundantes os calhaus resultantes 
da sua desagregação: ainda neste caso se observa 
concordância com a rocha encaixante e sub-vertica- 
lidade das camadas. 

Dentro de uma classificação morfológico-gené- 
tica, considera-se este afloramento scheelítico englo- 
bado no grupo dos jazigos estratóides de substituição 
metassomática em rochas sedimentares carbonatadas. 

A mineralização ocorre disseminada no seio da 
rocha escarnítica e mais raramente em finos leitos 
preferenciais, parecendo revelar tendência a associa- 
cões com idocrase. 

Numa perspectiva de valorização económica é 
provável estarmos em presença de ocorrências pobres 
não só pela irregularidade da mineralização como 
também pelo seu fino granulado. 
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Ill — PETROGRAFIA 


WI. 1. Rocha encaixante 


As bancadas de natureza calco-silicatada estão 
envolvidas por micaxisto por vezes profundamente 
alterado, O micaxisto revela textura granolepidoblás- 
tica e é constituído por quartzo, biotite e moscovite, 
como minerais essenciais, e acessoriamente por seri- 
cite, turmalina, apatite, plagioclase (Andesina An 38) 
e opacos, 

O quartzo é mineral frequente sob a forma de 
grãos com dimensão média de 0,2 mm, revelando 
fracturas normalmente ocupadas por óxidos de ferro 
que resultarão da descoloração muito profunda da 
biotite. Esta mica revela de facto uma perda acen- 
tuada de coloração passando a moscovite, mas em 
certas porções conserva ainda a cor original casta- 
nha avermelhada (ng=nm) variável até amarelo 
quase incolor (np). Inclusões de zircão com auréolas 
pleocroicas de tom castanho escuro. 

A moscovite parece poder ser atribuída a duas 
gerações: uma que teria resultado da biotite disposta 
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segundo a estratificação e outra posterior, cortando 
cristais biotíticos. 

A partir deste tipo xistento é possível estudar-se 
a passagem às bancadas calco-silicatadas já detectá- 
vel a nível de afloramento. 


Assim, desde o micaxisto onde a biotite é substi- 
tuída por clorite, evolui-se para uma banda essen- 
cialmente constituída por calcite com o aparecimento 
prévio de grãos de clinozoisite e de anfibola actino- 
lítica. A clinopiroxena diopaídica associa-se em peque- 
nos cristais (dimensão média de 0,87 mm), sendo 
incolor a muito levemente esverdeada, Este mineral 
com o quartzo e a plagioclase forma uma matriz 
granular que envolve esfena, apatite, grãos de miné- 
rios opacos e esporadicamente pequenas lamelas bio- 
títicas normalmente cloritizadas. 


A passagem às faixas calcíticas faz-se de modo 
brusco e é marcada pelo desenvolvimento de cristais 
de diopsido que são envolvidos por uma matriz resul- 
tante da profunda alteração da plagioclase onde se 
destacam grãos de clinozoisite, zoisite e calcite, assim 
como amígdalas de sílica criptocristalina. 


Il. 2. Rochas calco-silicatadas 


O aspecto macroscópico destes tipos litológicos 
mostra-se por vezes acentuadamente variável de aflo- 
ramento para afloramento e por vezes dentro da 
mesma camada, 


Em qualquer caso, no entanto, uma meteorização 
diferencial confere às rochas aspecto característico 
em que finas faixas mais salientes alternam com 
outras de aspecto reentrante. Este aspecto, associado 
à faixa exterior acastanhada de alteração constitui 
um factor auxiliar importante na prospecção de tais 
tipos litológicos. 

As rochas estudadas são de modo geral fina- 
mente granulares e bem foliadas, pelo que resultam 
bandas centimétricas de diferentes colorações com 
frequentes tonalidades esverdeadas alternantes com 
outras acastanhadas; no entanto algumas das rochas 
podem não mostrar litagem mas um aspecto maciço 
com tons de cinzento esverdeado. 


Uma mesma bancada pode revelar variações lo- 
cais com faixas essencialmente marcadas pelo desen- 
volvimento preferencial dum mineral como acontece 
com a idocrase, granada ou com minerais do grupo 
do epídoto. 

Desenvolvendo-se paralelamente à estratificação, 
veios pegmatíticos com espessura média de três cen- 
timetros podem evidenciar-se no seio das camadas 
calco-magnesianas, 


IH. 2.1. Estudo microscópico 


Os tipos litológicos de Rocamondo e F'reches 


mostram composição mineralógica idêntica pelo que 
passam a ser analisados em conjunto, 
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Sob o microscópio é visível a estrutura normal- 
mente bem foliada que deverá corresponder à estra- 
tificação da rocha original; a textura é granoblástica 
podendo passar a porfirogranoblástica, 

QUARTZO — Mineral comum a todas as rochas 
estudadas ocorre em grãos; eventualmente em bandas 
paralelas, Os cristais têm disposição em mosaico com 
os contornos de aspecto lobado e dimensões variáveis 
entre 0,03 e 2mm, Formas de contorno hexagonal 
surgem com carácter esporádico. 

PLAGIOCLASE — Cristais revelando tendência 
prismática com dimensões entre 0,04 e 2mm, São 
frequentes maclas da albite e associações albite-peri- 
clina, surgindo também geminações devidas a acções 
de natureza tectónica com terminação em bisel e 
planos de macla com fortes encurvamentos. 

Trata-se de Bytownite (Am 74). No entanto, al- 
guns dos tipos estudados revelam albitização mais 
ou menos profunda com desenvolvimento de Albite 
(An=8%). 

O mineral ocorre ainda intimamente associado ao 
quartzo formando uma estrutura em mosaico, cuja 
origem poderá estar relacionada com recristalizações 
após fragmentação periférica dos cristais mais desen- 
volvidos. 

Agregados de finas lamelas sericíticas associadas 
a grãos de clinozoizite e plagas calcíticas substituem 
com frequência os cristais plagioclásicos, podendo 
associar-se prenite que desenvolve uma rede de velos 
finos que cruzam os cristais. 

A. evolução da plagioclase no sentido de um fel- 
dspato potássico do tipo da microclina pôde ser obser- 
vada em algumas das rochas pelo desenvolvimento 
de manchas microclínicas nos cristais de plagioclase. 

CLINOPIROXENA — Sempre presente em todas 
as rochas sob a forma de grãos ou cristais com ten- 
dência prismática; as dimensões médias variam entre 
0,03 e 1,7 mm. 

São frequentes maclas polissintéticas muito finas 
com planos de macla encurvados. 

Sobre os cristais pode desenvolver-se uma anfi- 
bola actinolítica que mantêm a mesma continuidade 
óptica. 

Trata-se de uma piroxena da série do diópsido- 
-“hedenbergite com ângulo dos eixos ópticos (2V) va- 
riando entre 56º e 61º (valores mais frequentes da 
ordem dos 58º) e ângulos de extinção (c A ng) com- 
preendidos entre 42 e 45º, 

ANFIBOLA — Certos tipos calco-silicatados po- 
dem revelar o aparecimento duma anfibola actinoli- 
tica sob a forma de cristais primáticos, com termi- 
nações franjadas, dispostos alinhadamente com a fo- 
liação. Trata-se duma actinolite de cor verde leve- 
mente azulada (ng) variável até amarelo esverdeado 
(np) e com ângulos médios de extinção 17º (c A ng). 

IDOCRASE — Mineral que ocorre apenas nal- 
guns dos tipos estudados, mas onde toma grande 
desenvolvimento ocupando elevada percentagem em 
volume (Quadro I) (Fot, 1). 
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QUADRO | 


ANÁLISES MODAIS DE ALGUNS TIPOS CALCO-SILICATADOS 
(percentagens em volume) 


CcSi RO Cs6 KO Cs7 RO Cs1 RO CSA RO 

Quartzo 21,29 12,73 19,24 6,10 2,50 
Plagioclase 46,08 96,31 27,01 44,01 6,62 
Diopsido 26,59 24,51 28,50 42,25 9,92 
Actinolite 0,39 1,12 -— — — 
Idocrase — -— 22,02 = TS, 70 
Granada — — — e 0,56 
Esfena 2,19 2,82 1,39 2,36 0,17 
Clinoz + zoiz. 0,22 0,53 1,56 3,67 1,06 
Apatite (1,26 0,11 -— 0,11 0,11 
Prenite 1,93 0,56 — 1,05 0,56 | 
Úpacos 0,39 0,77 trs — mid 
Calcite 0,13 | Es = is = 
Fluorite — 0,29 — — E 
Scheelite E 0,25 0,28 0,45 trs 
Clorite — mm trs Fa er: 
Pontos contados 2279 | 2 832 2 888 2672 1794 


LOCALIZAÇÃO DAS AMOSTRAS 


CS1 RO—A S-40ºE e a 925 metros do marco geodésico de Lombo-Velho 


CS86 RO-- A S8-22ºE e a 1625 » %» > > > 
CS5S7 RO-—A S-22ºE e a 1650 » > ” » 
Cs10 RO—-A S- 4ºE e a 1175 > > > > 
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Fot. 1 — Cristal de idocrase com extinção em faixas atravessado por fino veio 
de quartzo. Luz pol. n+, (x25) 


As formas então reveladas assumem aspectos 
prismáticos ou quadrangulares com contornos rectilf- 
neos perfeitamente marcados pelo quartzo, Os cristais, 
com extinção característica por secções, atingem di- 
mensão máxima da ordem dos 4,4 mm, individuali- 
zando faixas bem definidas onde podem revelar as- 
pectos poincilíticos. 

GRANADA — Também a sua ocorrência está li- 
mitada apenas a alguns dos tipos em estudo, podendo 
aí ocupar grande percentagem sob a forma de cris- 
tais incolores corroendo os restantes componentes. 


Em todas as rochas estudadas e acessoriamente 
ocorrem esfena, zoizite, clinozoizite, (Fot. 2), apatite, 
opacos que em grande parte deverão ser pirrotite; 
com carácter mais esporádico calcite, biotite, clorite 
e prenite que pode surgir quer como produto de alte- 
ração da plagioclase, quer rodeando grãos de miné- 
rios, quer ainda formando veios que atravessam as 
rochas, 

Em quantidade variável, mas presente em quase 
todas as rochas a scheelite ocorre em pequenos grãos 


Fot. & — Aspecto textural de tipo calco-silicatado com cristais prismáticos de 
clinozoizite numa matriz quártzica. Luz pol. n+, (x62) 
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(só esporadicamente atingindo dimensões máximas da 
ordem dos 2,68 mm) de tom levemente acastanhado; 
surge inclusa na idocrase, ou pode ser envolvida pelo 
quartzo que a limita por vezes de modo rectilíneo 
pelo que é de admitir formação anterior a este mine- 
ral (Fots. 3 e 4), 


metamorfismo de contacto sob influência de tempe- 
raturas e pressões elevadas. 

Tais accões ter-se-iam exercido sobre primitivos 
leitos calcário-dolomíticos impuros intercalados na 
série xistosa, anteriormente afectados por acções me- 


tamórficas regionais, 


Fot. 3 — Cristal de scheelite no seio de piroxena e quartzo, Luz pol, n+, (x25) 


Fot. 4 — Cristal de scheelite intensamente fracturado envolvido por diópsido 
e quartzo. Luz pol., n+, (x25) 


IV — BREVES REFERÊNCIAS AOS ASPECTOS META- 
MÓRFICOS 


O desenvolvimento dos níveis calco-silicatados 


na vizinhança da rocha granitoíide aflorante na região 
leva a admitir origem relacionada com acções de 
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Sobre tal associação exerceram-se acções de me- 
tamorfismo de contacto impostas pelas massas gra- 
nitoides e primeiramente traduzidas por efeitos de 
carácter térmico; as associações minerais presentes 
nos tipos litológicos estudados e particularmente a 
ausência de volastonite fazem admitir influência de 
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temperaturas variáveis entre 500 e 600% (Zharikov, 
V. A, 1969). 

Em parte com carácter simultâneo e posterior- 
mente âqueles efeitos desenvolveram-se acções de me- 
tassomatismo com acarreio de substâncias. Neste 
caso foram importantes os fluidos responsáveis pelo 
transporte de Si, W e Fe, além de fluor, enxofre a 
fósforo, estes condicionando o desenvolvimento de 
fluorite, pirrotite e apatite; também o surgimento da 
vesuvianite poderá estar relacionado com o acarreio 
de fluor, embora a sua formação deva estar condicio- 
nada por acções metamórficas de contacto. 

Fenómenos regressivos ligados a acções de ca- 
rácter hidrotermal tiveram também influência na 
evolução dos tipos calco-silicatados, como pode ser 
atestado pela substituição das plagloclases por seri- 
cite e clinozoizite. 

De acordo com os conceitos de fácies de alcali- 
nidade (Zharikov, V. A., 1969) as rochas estudadas 
serão incluídas na fácies do tipo piroxena-plagioclase, 
enquanto que atendendo à composição mineralógica 
se incluirão na fácies das corneanas piroxénicas ('Tur- 
ner, F. e Verhoogen, J. 1951) sendo as seguintes as 
principais associações mineralógicas: 


quartzo + plagioclase + diopsido + clinozoizite 
quartzo + plagioclase + diopsido + actinolite + cli- 
nozoizite 


quartzo + plagioclase + idocrase + granada + clino- 
zoizite 
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Método directo de distribuição de momentos para o cálculo 
de estruturas sujeitas a esforços horizontais 


RESUMO 


Neste trabalho deduzem-se fórmulas que aplica- 
das ao método de distribuição de momentos em 
consola reduzem este a um ciclo para cada nó da 
estrutura fazendo desaparecer o inconveniente da 
lenta convergência. 


INTRODUÇÃO 


Em publicação anterior [1] abordou-se o proble- 
ma das estruturas sujeitas a esforços horizontais nos 
seguintes casos: 


1. Pórticos simétricos 
2. Pórticos obedecendo ao princípio dos múltiplos 
3. Pórticos assimétricos 


Em [1] utilizou-se o método de Grinter — «Canti- 

lever Moment Distribution: [2] que consiste em mo- 
K 

dificar a rigidez dos pilares para : e o coeficiente 

de transmissão de momentos para —1.0, Isto justifi- 

ca-se visto a rotação numa extremidade duma con- 


sola sob a acção dum momento M ser e o coefi- 


ciente de transmissão ser —1.0 para a convenção 
ML 
4EI 
numa viga com apoio simples sob a acção do mesmo 
momento. 

Em [2] aplica-se o mesmo processo depois da 
decomposição da estrutura em pórticos simples simé- 
tricos. 

Em [3] substituiu-se a estrutura por uma consola 
equivalente e aplicou-se o método de Grinter a esta 
estrutura de maneira a obter valores aproximados 
das rotações e dos ângulos de deslocamento nos nós. 
A partir daqui com as equações de deformação e 
por iteração muito rápida obtém-se os valores finais 
que permitem determinar os momentos na estrutura. 


usual de sinais, comparada com a rotação 


Manuscrito recebido para publicação em 28/1/7171. 
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SUMMARY 


The paper works a few formulae to reduce the 
cantilever moment distribution methode to a cycle 
for each joint of the structure with sidesway and 
overcoming the inconvenience of the sloro conver- 
gence. 


Ora quando a rigidez das vigas é pequena em 
relação à rigidez dos pilares, caso muito frequente 
em edifícios de vários pisos, a convergência no Mé- 
todo de Grinter é muito lenta e, por conseguinte, 
vão-se determinar fórmulas que permitem para cada 
momento actuante num nó depois de fazer o equilí- 
brio transmitir momentos finais aos outros nós da 
estrutura. 


MÉTODO DIRECTO PARA NÓS FIXOS 


Consideremos uma barra ik, numa estrutura em 
que os nós não se deslocam 


fig. 1 


e fazemos actuar na extremidade i o momento M ne 
Escrevendo as equações de deformação 
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(1) 
M, = Ki (op TT Pi) 
em que 
4 EI 
K, = ———+ (2) 
l 
ix 
é a rigidez absoluta da barra lk. 
Mas M, = Ci Pi 
(3) 
— Mo = Cu Pk 


sendo C, e OC, os factores de rigidez da barra ik 
respectivamente nos nós ie k. 
Substituindo q, de (3) em (1) 


Dae M, 
M,, = K,, (e; — 2€ ) (4) 
ki 
M,. P 
ME, (———— + —) 
ki 1k C. 2 


Eliminando M,, em (4) obtém-se 


0,75 K,. +C. 

M,JDK,,—-— q (5) 
ik ik Pi 
K,, + Ci 
e então 

0,75 EK, FC 
Pe gg es do ; E (6) 

K me Cu 


Eliminando y, em (4) obtém-se 


1 
A re e — M,,=Ty My (7) 
2 0Q5KR, + €C. 
1 Cc | 
e = — e (8) 
2 BMb Rj + Cu 
é o coeficiente de transmissão de momentos. 
A fórmula (6) pode-se escrever 
Cc, CE Vik K (9) 


sendo v,. o coeficiente de restrição da barra ik na 
sua extremidade k 


Ci 
0,75 + — 
ki (10) 
FpD e 
Ko 
1 
K 


ik 


e varia de 0,75 (extremidade simplesmente apoiada) 
a 1 (extremidade encastrada). 


O coeficiente de transmissão r,, em função de v,. 
será 
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1,5 


v 


Pi 2—— (11) 


ik 


As fórmulas 6, 7, 8, 9, 10 e 11 permitem resolver 
uma estrutura com nós não deslocáveis apenas com 
um ciclo. 


MÉTODO DIRECTO PARA NÓS DESLOCÁVEIS 


Consideremos agora uma barra AB com um 
momento M, actuando em B em que os nós se des- 
locam 


M=C or (12) 
M=C Pp (13) 
Cha CpotCay (36) 
C,y é O factor de rigidez da viga. 
Mas E Sao sam 
e 6 pelo teorema de Mohr é igual a 
M lp +4M 4C, nPa 
p= É me (16) 
AR Kan Kan 
substituindo este valor em (15) 
CanPa  4CsanPa (17) 
à À Ma a 
Coa Kan 
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Dividindo por q, 


, Ca - 4 Cb Ea 
Coa K ap 
K 
AB 
Fazendo Ko = (19) 
AB 4 


e substituindo em (18) obtém-se 


SK ———— (20) 


A fórmula (20) dã o factor de rigidez numa 
extremidade em função do factor de rigidez da outra 
extremidade e da rigidez modificada que juntamente 
com o coeficiente de transmissão — 1 permite resol- 
ver o pórtico simétrico pelo método directo. 

Vai-se exemplificar primeiro pelo método de 
Grinter e depois pelo método directo para o pórtico 
da figura 3. 


Ê a 
E K=20 
K = 1.6 12 
6 T 
eg y pigs: 
D* 
K = 20 
| K =2 10 
5 T 
RE = so 
E K = 2.4 
K=3 g 
B' 
K=1.2 
K=3 8 
A a 
|| 
Fig. 3 


METODO DE GRINTER 
Devido à antimetria a rigidez das vigas é 1,6 Kv. 


Assim a rigidez modificada encontra-se representada 
na fig. 4. 


MOMENTOS FINAIS 


iza 804 fato 450 + 061-—007|+0.02 0.00 /|+ 34.31 
K' = 0.4 
+30 — 825/+ 45 — 061 |+0.07— 0.02 | + 25.69 
di | | dó ibid dididi + 013 — 0.08 || + 40.88 
K = 05 
+40 — 10.42] + 10.81 — 1,43 |+ 0.76 — 0.13 |+ 0.03 — 0.01 ||+ 39.11 
pia 66 pi di | tai pr | po a 
K' = 0.75 
+ 42 —19.09/+ 1563 — 3,55 |+ 2.140,48 | + 0.19 — 0,04 + 36.80 
B 
ii +42 — 19.09 — 3.55 — 048 — 0.04) + 18.84 
K' = 0.75 
| + 42 — 19,09 + 3.55 + 0.48 + 0.04] + 65.16 
m + 
Fig. 4 
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MÉTODO DIRECTO a PO TUDO: 


C =04 
ED 0.4 + 3.45 
Aplicando a fórmula (20) obtém-se 
3.0 | 
Do Co =04-—-—————— = 0,35 
Cc =05————— =0,75 DE 0.4 + 8.0 
BA 0,75 + 
3.0 + 0.85 
CB 07,5 +2.55 - didi 
; Ê 36 + 0,44 
Cc =065 SEU = 0,45 Cc =975 o. = 0.63 
CD 0.5 + 4.18 BC 0.75 + 4.04 
o E No D No' € No B MOMENTOS FINAIS 
| 
0.36  +30— 3.21 + 6.45 + 082. + 0.19 | + 34.25 
0.35 + 3.21 | + 30 — 6.45 — 0.82 — 0.19 || + 25.75 
3.0 | 
== | 
0.45 | — 0.46 | +40— 8.29. 47.81 + 1.81 || + 40.87 
| | | 
0.44 + 0,46 + 8.29 | +40 — T.81 — 1.81 | + 39.13 
3.6 
E 
0.58 — 0.06 — 1.15 | 4 42 — 10.29 + 16.64 || + 4714 
| 
| | 
0.63 + 0.06 +115. +10.29 | +42 1664 | + 36.86 
1.8 | 
A "B | 
0.75 — 0,02 — 0.34 — 3.03 | + 42 19.81 || + 18.80 
E + 0,02 + 10.34 + 3.03 + 42 + 19.81 || + 65.20 
Tr * | | 
Fig. 5 
| 0.36 nom F | , 
No E Equilíbrio — — x 30-32] dão os rRorentos de encastramento para os deslo 
3.36 camentos reais 4; na 2." faz a distribuição destes 
momentos utilizando o método directo para o caso 
o 045 de nós não deslocáveis, 
M pe = -— 3.21 X “345 — — 0,46 
REFERÊNCIAS 
M...=— 0.46 X - E 0.06 
CB | 418 | [1] Cálculo de Estruturas sujeitas a esforços hori- 
zontais. A, de Pinho Morgado, BOLETIM DA 
— 07.5 SOCIEDADE DE ESTUDOS DE MOÇAMBIQUE 
M.,-— 0.06 xXx -——— = — 0. 4 
BA X 255 ua — N.º 120 — 1960. 
[2] Joit translation by cantilever moment distribu- 
Para os outros nós procede-se da mesma ma- tion. L. E. Grinter and €. H. Tsad, TRANSAC- 
neira e o desenvolvimento no cálculo apresenta-se TIONS AMERICAN SOCIETY OF CIVIL EN- 
na fig. 5. GINEERING, VOL. 118 — 1954, 
RYGOL resolve este mesmo exemplo em 2 fases: [3j Structural Analysis by direct moment distribu- 
na 1.º determina os factores de deslocamento que tion, 3. Rygol. 
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Betão normal reforçado com armadura metálica 


7 — «ECONOMIA — RENDIMENTO DUMA SECÇÃO» 
7.1 — Generalidades 


Face a um problema concreto importa definir 
a secção mais económica ou seja, em principio, a 
que exige a menor quantidade de materiais. 

No que se refere à armadura esta está directa- 
mente relacionada com a altura da secção que, por 
sua vez, pode ser condicionada ou por condições cons- 
trutivas ou por condições de deformabilidade. 


71.2 — Altura da secção (h,) 


Dão-se seguidamente algumas normas que per- 
mitem determinar a altura total aproximada da sec- 
ção, de modo a permitir estabelecer rapidamente, 
em «ante-projecto», este elemento geométrico, em 
função do vão 1. Para mais detalhes ver também 
capítulo em que se estuda a «Deformação». 

Em 7.21. e 7.2.2. estudam-se casos correntes em 
que a > 0,5; em 7.2.3 casos em que q < 0,5. 


7.2.1 — Caso de pontes 


a) Tramo simples (com I< 20m) 


— lajes 
o <h < (7.1) 
— vigas (com I > 20m) 
E o qu (7.2) 
20 20 
b) Tramos contínuos 
— lajes 
h, > = (7.3) 
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— vigas 


(7.4) 


7.2.2 — Casos correntes 


Em casos correntes de estruturas em geral, pode 
tomar-se genericamente para fins de «ante-projecto» 


(7.5) 


7.2.3 — Regulamento nacional 


As indicações dos Art.o 53 e Art.” 63 do R.E.B.A. 
fornecem metodologia utilizável para a determinação 
da altura útil das secções rectangulares e, corres- 
pondentemente, de h, atendendo à deformação. 


7.3 — Espessura da alma (bo) 


Seguem-se alguns valores para o caso de q > 0,5. 
O valor médio corrente em pontes é da ordem de: 


, = 20 em (7.6) 
O valor mínimo no centro do vão: 
h 
bo = 9 cm + Delem) (7.7) 
40 


Quando se adapta este valor mínimo, haverá que 
aumentar progressivamente a espessura no último 
quarto do vão de modo que se obtenha nos apoios 
uma espessura suficiente para resistir aos esforços 
de corte e às tensões locais de contacto com as 
ancoragens de pré-esforço, quando existam. 
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7.4— A secção económica. Localização do Problema 


Dado o valor do momento máximo M, e fixado h,, 
obtém-se para ordem de grandeza da força total de 
tracção, o valor: 


F= (7.8) 


A =n bh (7.9) 


e, de (5.50) 


c=(,)—(p.6) 


De (5.11), é: 
| : bh 
ato 2 , bg * - Voa 
Ou 
| ' - am Ma 
e 4 td [ad 2“ - pes 
2 bs oz | Ea h 
Ê 
fp M, 
ey - >< 
É - " U,5 “ 5a “a h 
donde 
n M 
E qt. maça! (7.10) 
ha by f h : 


sendo, conforme (5.40) 


1 uh = 0,5 «E é 4 E (7.11) 


z 


Assim do ponto de vista da armadura, 
— ...à economia depende directamente de 


E «h 


Do ponto de vista do betão, 
— ...& economia depende da grandeza 


M, 


dp = = | ; (7.12) 


Ê 
e de ua 


7.5 — Análise referente ao betão 
Quanto menor for 3. e maior for o',, menor será 


a área necessária de betão A, e, portanto, mais 
económica do ponto de vista do betão. 


2 


A 5, chamaremos «coeficiente de rendimento» 


da secção. 
Por outro lado (7.10) permite obter b. 


3 M 
h E o r Ea] a Aa 
! " | = hi (1.13) 
b ha 
5 
A (=) chamaremos «coeficiente de flexão» 
b 
da secção. 


Com efeito numa secção rectangular, de material 
homogéneo, este coeficiente tem o valor 6. 


M. —M, 
Sa qe Em 
donde 
M 
b = 6 E eae ny 
ps o hº 


De (7.10) pode obter-se, A,, mais explicitamente: 


A =. (7.14) 


7.5.1 — Coeficientes de “rendimento” da secção 8, 


E] ô, hd 
e de “flexão” ( ) 
Tp 


Determinam-se seguidamente vários valores de 8, 
função da geometria da secção, n, e m,, da posição 
da linha neutra q e de diferentes valores do coefi- 
ciente (q. €,). 


E: 
ri 
5, ai b am EM va RD. 
a 0,5xXax0,,| (8,)—(g.8) | 
(7.45 
para 


n, — Ver (5.3); 

(E,) — ver (6.28) e (5.29); 
(p-e£,)) —ver (5.23) e (5.30); 
Oo; Ver (5.12); 

wa — posição de 1.n. 


Quanto menor for «5,» mais económica será a 


secção do ponto de vista do betão. 


dr 
o 


Os quadros 7.1 e 7.2 dão valores de ( 


) e de 


(8,), respectivamente para (p.e )=0 e e = (e), 
e (gy ED)= 0,10 e E (e )— 0,10 
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0,40] 1,39 


Tal 


——e À) —eeee a À e | ——— À 


5,55] 4,46] 3,93] 11,44 
1,91] 1,91] 201] 1,14 


O À =D TT O TT E O a ana 
- > = O e 
—— — dd —— os] o 


5,18] 4,13 
1,77] 1,77 


3,99| 10,93 
1,84! 1,09 


ô 
Fm — em cima; “8,” — em baixo. Os valores (7) são os inversos de m, 
b 


QUADRO 7.1 
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0,585 | 0,635 


0,343 | 0,418 | 0,493 


13,53 | 12,05 
7,91 | 7,65 


14,52 | 12,05 
4,98 | 5,04 


10,42 8,76 
6,10 5,06 


8,79 7,19 
o,1d 4,97 


5,56 | 10,51 
3,16 | 2,75 | 


5,02 | 10,05 
285 | 263 


6,02 5,03 
3,52 3,20 


4,43 | 9,26 
252 | 2,43 


fi E E. 


7,57 | 5,58 | 
260 | 233 


4,05 | 8,47 
230 | 2,22 


ô, E pá 
(=) — em cima, à —em baixo. Os valores (O são os inversos de m, 


QUADRO 7.2 
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Verifica-se que 5, aumenta com a diminuição de a 


e com o aumento de (y .Ee,); é maior para as secções 
rectangulares que para secções T. 

Valores de n, pequenos não são normalmente os 
mais económicos, 

Para n, da ordem de 0,20 a 0,30 obtém-se valores 
de 3, pequenos e perto dos mínimos. 


Pe at: 


Para cada a os coeficientes ( 


b 
variação sensível para m, < 0,30 e n, > 0,20. 


+ E 3 
No Quadro 7.3 apresenta-se uma síntese do 
ô 
exposto, indicando ordens de grandeza de (—L) de 
% 


modo a determinar aproximadamente valores de b 


que conduzem à secção económica do ponto de vista 
de betão. 


«SECÇÃO ECONÓMICA» — Definição da largura b — 
(8,/m,) 


m = 15 


8,3 
B. ne o 
0,40 58 | 65 
48 5,5 
0,80 4,2 Fr 
e A 


3,4 4,0 


3,0 


Os valores deste quadro conduzem à definição 
aproximada da largura b da secção mais económica 
do ponto de vista de betão (A, — mínimo). 


e E o e em a 


QUADRO 7.3 


u econ. 


Para pré-dimensionar uma secção do ponto de 
vista do betão procede-se da seguinte forma: 


— patima-se M,; fixa-se h, e Che? 


— Pelo Quadro 7.3 retira-se (8,/n,) de de (7.13) 
obtém-se a largura b; 


—b, será o menor possível condicionado por 
condições construtivas e pelo esforço trans- 
verso; 


—para h, tomar-se-á um valor da ordem 
de 0,20 a 0,30 de h. 


7.6 — Valores de «u» económicos 


A finalizar este capítulo e como elemento infor- 
mativo, apresentam-se seguidamente ordens de gran- 
deza dos valores de q correspondentes a secções eco- 


nómicas do ponto de vista do betão (A, — mínimo). 


É, como se sabe, u= 
b.h? 


de (7.13), obtém-se: 


(7.16) 


Assim, tomando para (5,/n,) os valores médios 
do Quadro 7.3, obtém-se os valores de do Qua- 
dro 7.4. Note-se que, por questões de segurança à 
rotura, valores de « maiores que os indicados para 


«=0,80 serão dificilmente atingíveis. 


o',,— tensão admissível num dimensionamento por 


«tensões de segurança». 


115 | 130 | 150 | 170 | 200 | 60 | 


QUADRO 7.4 
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T a 
200 
43 
40 | 46 
M, 
"ph 
55 


7.7 — Relação (l/h,) económica 


O estabelecimento das relações 1/h, mais cor- 
rentes e para as quais já se indicaram alguns crité- 
rios em T.2.1, 7.2.2 e 7.23 importa ser visto em mais 
detalhe. 

Esta relação está dependente dos condicionalis- 
mos próprios da obra em estudo, da deformabilidade 
aceitável para a obra, da procura de soluções econó- 
micas e da largura b ou b, necessárias. Estes últi- 


mos elementos (b ou b,) estão por sua vez relacio- 


nados com problemas de esforço transverso e de 
deformabilidade transversal da peça em estudo. 

É fácil demonstrar que tanto a quantidade de 
armadura A como a área de betão A, necessárias, 


são tanto menores quanto menor for a relação (1/h,). 


Com efeito: 
M | 
à sra Eae af po tado 
gh Ea Tum h, h 
io 
j M 
A,=a,.bh E qe 
(Pa h 
1 , | 
mo coa e 3 tt À (718) 
“bs Pia h, h 


Simplesmente, de (7.13) verifica-st que aumen- 
tando h, ou h haverá que reduzir b para se obter a 
a secção económica e este valor b não pode descer 
abaixo de certos valores por questões de execução, 
de esforço transverso e de deformabilidade trans- 
versal, como se referiu. 


Assim, fazendo; 


7.19) 


D. M, 
kp — Ny ha ne > 


O momento total M, pode obter-se, como é 
sabido, a partir de expressões do tipo: 


º q. 
M, = Em (7.20) 
donde: 
| - (— + A su ( h, ) 
h, o LU “mm Fa, , h / 
56 


e para, h,/n=1,07 


| É à nincê / “ui “AM l 
———= | == Udo) es o Ag 
oh, V KR, 4% f%) 4 


Il —vão da viga; 
q— carga total por unidade de comprimento; 
h, —altura total da viga. 


A expressão (7.21) é geral e pode ser aplicada 
a qualquer caso concreto. 


hoy 


D) 


Os vatores ( E ) são tirados do quadro 7.1 ou 


h, 
do quadro T.2. 


»od 


Alguns regulamentos relacionam a tensão admis- 


sível O ha com o valor de k,. 


Estão neste caso as “normas inglesas”. 

As “normas americanas” indicam que, para vigas 
isoladas, k, deve ser menor ou igual a 32, 

Pode dizer-se que, para vigas isoladas que sejam 
solidarizadas entre si, é corrente tomar-se valores 
de k, da ordem de 40 (k =40). 

Sobre tudo em vigas pré-fabricadas, atendendo 
ao seu transporte e montagem, não convirá utilizar 
valores de k, apreciavelmente maiores que o valor 


indicado. Ver também 7.11. 
+ do mk 
Em vigas isostáticas é, como é bem sabido, —8 


e a expressão (7.21), para k, =40, pode então escre- 
ver-se: 


(7.22) 


No Quadro 7.5 e no Quadro 7.6 apresentam-se 


t 


ordens de grandeza da relação 5) 


— O Quadro 7.5 refere-se a secções de betão 
armado corrente (B. A.) tendo-se tomado: 


«=0,35/0,40 : o',.=80 kgf. cm-2 


1) 
para secções E — ( - ) =6,0; 


LM 


Ê 
para secções T — ( É =T,, 


o 
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— O Quadro 7.6 refere-se a secção de betão 
armado/pré-esforçado e betão pré-esforçado, (B.A./P.) 
e (B. P.), tendo-se tomado: 


04=0,80 o',,— 120 kgf. ecm-2 


j 


3 
para secção E — ! . ) = 3,4 


yo 


Ar 
para secção T — ( - ) =4,5 
"it 


4=0,35/040 
o',.=80 kgf.cm-2 


| 
— 40 
b 


QUADRO T.6 
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Como se verifica os valores do Quadro 7.6 estão 
concordantes com os valores indicados em T.2.1 e 
7.2.2. 


Note-se que se tomou 1 = 40; para outros valo- 


res desta relação obter-se-iam valores um pouco 
diferentes. 


** + 


Os valores obtidos no Quadro 7.5 e no Quadro 7.6 
serão comparados com os valores da mesma rela- 


TR 

ção 5) impostos por questões de deformação e 
L 

que serão estudados em capítulo próprio. 


71.8 — Considerações quanto à rotura 


As secções devem ser também verificadas à 
rotura. Para w > 0,8 é este dimensionamento que, em 
geral, condiciona a secção. 

Este estudo será feito em capítulo próprio. 

Refere-se, no entanto, que para anteprojectos 
se podem, em geral manter os critérios apontados 
neste Capítulo 7 desde que, no que se refere à secção 
de betão, se tomem valores não inferiores aos indi- 
cados para q=0,8. 

Com este critério fica normalmente, garantida 
a segurança à rotura, 


7.9 — Análise referente à armadura 


De (7.8) verifica-se que a quantidade de arma- 
dura está directamente relacionada com 


(7.23) 


sendo E =(E )-(q.E,) 


M, — é o momento actuante; 
h— é a altura útil como definida na Fig. 5.1. 


Como valores aproximados de “e” indicam-se 
[ud 


os do Quadro 77 na secção mais esforçada. 


“e nr 
ad 
B.A. B.A./P. B.P 
02205 | 06a 08 1.07 
0,83 0,75 0,60 
0,86 080 | oq 
QUADRO 77 


Este quadro e a expressão (7.23) permitem o 
predimensionamento imediato da força de tracção a 
mobilizar na armadura. 


OT 


710 — Exemplos 


7.10.1 — Dimensionar para fins de ante-projecto a 


secção central de uma viga de betão armado 
(B. A) com um vão |=20 m, submetida a uma 
carga uniformemente distribuída q=3,0.t/m 


Admite-se o”, = 130 kgf cm-? e a, = 2400 kgf. cm-2. 


De (7.13) e (7.14), obtém-se: 


) M 
AM == = e - 
“os h 
É: 
M 3 >< 20º = 150 tf.m 
' ' 
Toma-se: 
h= 1 =20/20=10m e b=9cm 
2 


Em B.A. é (p.e)=0 ee =(e ) 
| Led ti 


a) — Secção E 


Do Quadro 7.3 toma-se ( Ê ) =5,8; 
= od 
será: 
b = 5,8 x 150108 73 em 
130 x 06º 
será: 


A,=b ; h, =73x100=7300 em: 


Do Quadro 7.1 é “8,” = 6,10, obtendo-se: 


, 610 150><105 
Cb 130 ** 96 si 


Este valor é muito semelhante ao valor real 
tomado A = 7300 cms. 


O valor de qu será: 


valor que se integra nas ordens de grandeza dos 


valores indicados no Quadro 7.4 para q entre 0,40 
e 0,50. 


o8 


— No que se refere à armadura, obtém-se, Qua- 
dro 7.7: 


M 50 se 108 
po to 10x<I6 oo asdugf 


Eh 0,83 x 96 


A = 188.250/2.400 = 78,0 em 


4 + x 


— Num dimensionamento directo par tensões de 
segurança obtinha-se: 


A =1,10 %x73xX96=77,1 em? 


oa 


-— Num dimensionamento à rotura (ver tabelas 
do Autor), obtinha-se: 


n*=1,5X22,3=33,5 


A =1,08 %xX73xX96=75,7 cm 


b) — Secção T 


Do Quadro 7.3 e de (7.13)), tem-se, para «= 0,40: 


ae 15050108 


b = 6, 
130 x 96º 


Fazendo h,=0,20X96= — 20 em 


e b,=0,25x78=20 cm 
A, = 3169 em? 
Do Quadro 7.1, obtinha-se: 


2.70 >< 150 >< 105 | 
A, = DDD = 345cm 
130 >< 06 


— Para armadura, do Quadro 7.7, tinha-se 


150 >< 105 
PF. = me SL Oo RE f 
nt 0,86><96 cida 


A. = Fat: 2400 = 757 em 


— Num dimensionamento directo por tensões de 


segurança obtinha-se: 


vo 150105 1 2087 
o 7896 * 2400 2400 000% 
o. = 0,0087 = QUI 
A f a 0,87 , 
A, = 0,01 x 18>(06 = 74,9 cm? 
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o. 


— Num dimensionamento à rotura (tabelas do 
Autor), obtinha-se: 


à — 15x 150x 106 


= 0,157 
78 >< 96º x 200 


1,5 >< 150 x 105 


fl EE sofia Ss aque Fani 
* 0,89 96» 3.600 


— 732 em? 


7.10.2 — O mesmo problema em betão pré-esforçado 
total (7,,—0), B. P. 


150 = 105 


ces Ei 
E AR Ta e 06 


= 62,0 cm 


Toma-se b=64 cm; b =20 cm; h,=20 em 


150 >< 105 
qi e 22: DHASÔ US 
“ = 0,74 96 tuas 


Fat = 211 tf 
Adoptando a secção indicada na Fig. 7.1, tem-se: 


A, = 2.880 cm: 
1=2,67xX 108 em? 


A 211.000 + 56><2 11.000 318 — 
4 2880 — 267><108 A L 
— 15) pd 105 ITS = — 118,3 kg Ei cm 
ME = + - 
+ 2,67 x 108 a PO 7 =+09 » 


Os valores obtidos são concordantes com os 
valores impostos, atendendo a que se tomaram os 
valores correspondentes a «=0,8 e não q=1,07 por 
razões de segurança à rotura. 

A armadura necessária, admitindo O = 10 000 


kgf. cm-?, será da ordem de: 


A = 211.000/10.000 = 21,1 cm? 
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710.3 — O mesmo problema em betão armado pré- 
-esforçado B. A.JP. 


Tomando « —1000 kgf. em-? e kap=1,0, e usando 
a geometria da Fig. 7.1, obtém-se: 


ve à .u 

[ (10.009 + 1,000) + 1.000] q, = ——— 
[e 

à = 0,287 % e 


À = A. = 16,9 cm: 


Para comparação apresentam-se resumidamente 
os valores obtidos: 


—de 7.10.1.0) — B. A. (Sec. [)—7.500 em; 76 cm? 
—de 7. 10.1.b) — B.A. (Sec. 7)—3.200 cm; 74 cm? 
—de 7.102 —B.P. (Sec. T)-2.880 cm”; 22 cmº 
—de 740.3) —B.A./P. (Sec. T)—2.880 em?; 17417 cm? 


710.4 — Determinar a ordem de grandeza do maior 
vão possivel de uma viga reclangular de 
betão armado simplesmente apoiada, subme- 
tida ao seu peso próprio e a uma sobrecarga 
equivalente a 50% do seu peso próprio. 


0',.=89) kgf. em-2 


Tomando 
[ M. 
Ay = 5,8 >< > 
bi bs h 
2500= A, ] 
M = 6x ————— se — se | 5x h 
= 2X 5210010000 *"h, 
5.8 2500 x A, 
A =— 15 = ———— se 15 >< 1x 1,07 
A o OOo 2 
l== 18m 


7.10.5 — Resolver o mesmo problema para uma viga T 
em betão pré-esforçado total (7,, = 0) 


o',,=140 kgf. cm-2 


Tomando 


Ss 


T bs h 


A b —— 2,10 De 


2500 x A, | 
E elo h 


fas dE SE jm 
” * 100><8>< 10.000 *'h, 


9 


2,10 2.900 = Ar, 
a erre SS e ODE 
140 * 8 e 1000 x< 10.000 | 


1=67m 


Os valores obtidos de | são valores teóricos, 
Para uma secção rectangular obtinham-se valores 
de 1 da ordem de 30m. 


7.11 — Considerações Genéricas sobre “k," 


Nos exemplos anteriores considerou-se k —40 
como indicação de ordens de grandeza e de valores 
médios e não, naturalmente, como valor único, 


60 


Com efeito, em certos casos, como por exemplo 
vigas de cobertura, podem atingir-se valores de K, 
da ordem de 70 ou mesmo maiores, 

Noutros casos, como por exemplo em pontes, 
pode interessar aumentar a largura b para cons- 
tituir, no todo ou em parte, o próprio tabuleiro. Nes- 
tes casos k, toma, por vezes, valores da ordem de 


20, ou menores. 


Para valores de k, =70 os valores da relação 


t 
serão obtidos dos indicados nos quadros dividindo 
aqueles valores por 1,32; para k, — 23 bastará multi- 
plicar aqueles valores por 1,32 para obter cs corres- 


pondentes valores de 
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Resumo dos Artigos publicados na «Técnica» n.º 443 


— 


C. D. U. 378.962.4:550.8[624, 628] 


L. AIRES-BARROS 
3. F. PINHARANDA 


O ENSINO DA GEOLOGIA NO |. S. T. PARA O CURSO 
DE ENGENHARIA CIVIL, EVOLUÇÃO E VISÃO PROSPECTIVA 


Técnica N.º 443 — LI|— 101977, p. 1a 6 


Os autores apresentam um resumo da evolução do ensino 
da Geologia para estudantes do Curso de Engenharia Civil 
do |. 8. T. nos últimos 20 anos. 

Em seguida é apresentada uma visão prospectiva do que 
deverá ser esse ensino e do equilibrio a haver entre os 
conhecimentos dos fundamentos da Geologia, das suas me- 
todologias e das suas aplicações a problemas concretos da 
Engenharia Civil. 

às principais conclusões são apresentadas sob a forma 
de um diagrama onde se doseia a influência e a sobreposição 
que hã entre os conhecimentos importantes nas diferentes 
fases da resolução de problemas reais da Engenharia Civil 
fortemente dependentes da Geologia. 

São feitas, ainda, algumas considerações sobre os tipos 
de exames que os alunos realizam. 


a = a 


C. D. U. 627.62:627.89 
J. LAGINHA SERAFIM 
DESASTRES DE BARRAGENS EM 1976. O CASO DE TETON 
Técnica N.º 443 — LI] — 10-1977, p. 7 à 17 


Este trabalho trata da rotura de barragens durante o ano 
de 1976 e está na sequência do trabalho apresentado pelo 
autor sobre o mesmo assunto no Laboratório Nacional de 
Engenharia Civil em 1960. Apresentam-se descrições de rotu- 
ras em barragens relativamente pequenas (mais de 15m de 
altura) na Carolina do Norte, no Paquistão e na Noruega, 
mas o trabalho trata essencialmente da rotura da barragem 
de Teton nos Estados Unidos, em Junho de 1976. Esta é a 
barragem mais alta (103 m) que até hoje rompeu. Apresenta-se 
a sequência dos acontecimentos durante a rotura, atribuindo 
esta ao facto da rocha ser muito permeável, à Ineficiência 
das cortinas de injecção, e ainda ao facto de existir uma 
trincheira corta-águas, na qual se produziu um efeito arco 
que possivelmente permitiu fracturação hidráulica. À erosão 
interna e subsequente erosão no aterro perto da rocha foram, 
sem dúvida, as causas principais. 
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Ano Lil — Outubro 1977 


C. D. U. 621.380:537.311 


H. J. 8. ABREU SANTOS 


SIMULAÇÃO NUMERICA DAS ESTATÍSTICAS DE MAXWELL- 
-BOLTZMANN, BOSE-EINSTEIN E FERMI-DIRAG, BASEADA 
NUM GERADOR DE NUMEROS PSEUDOALEATÓRIOS [1] 


Técnica N.º 443 — LIL — 101977, p. 19 a 38 


Faz-se um estudo comparativo da obtenção por simulação 
numérica, no dominio das pequenas populações de particulas, 
das funções de distribuição correspondentes às estatisticas 
de Maxwell-Boltzmann, Bose-Einsteln e Fermi-Dirac. É uma 
aplicação dos métodos de Monte-Carlo à simulação de um 
processo não determinístico, o das colisões aleatórias de 
uma população de particulas, utilizando um gerador de nú- 
meros pseudoaleatórios. Estuda-se a obtenção do equilibrio 
termodinâmico para populações correspondentes a uma só 
estatística, do tipo da lei de radiação de Planck, e também 
o interacção entre particulas de estatisticas diferentes, como 
electrões e fonões. 

Este trabalho é um prolongamento de [1]. 


ET o ET a E 


C D. U. 552.1.067.1 


A. F. FERREIRA PINTO 
A. F. BARROS 


OCORRÊNCIA DE SCHEELITE EM ROCHAS CALCO-SILI- 
CATADAS DA REGIÃO DE VILA FRANCA DAS NAVES 


Técnica N.º 443 — Lil — 101977, p. 39 a 45 


A presente nota refere a ocorrência de scheelite em 
bancadas de rochas calco-silicatadas das zonas de Roca- 
mondo e de Freches (carta geológica 18-A). Estas rochas 
constituem um conjunto grosseiramente paralelo com direc- 
ções N70ºW verticais a sub-verticais, 

Mineralogicamente trata-se de rochas com textura grano- 
blástica a porfiroblástica essencialmente constituídas por pla- 
gioclase, clinopiroxena diopsidica e anfibola, a que se asso- 
ciam idocrese e granada. Acessoriamente estão presentes 
esfena, zoizite, clinozolzite, apatite e opacos. A scheelite 
ocorre em quase todas as rochas em pequenos grãos. 

à formação e evolução das rochas calco-silicatadas estará 
relacionada com acções de. metamorfismo de contacto exer- 
cidas pelo granito no-szio do qual se localizam os retalhos 
xistosos que englobam os niveis estudados. 

A mineralização scheelítica ocorre normalmente dissemi- 


. nada no seio da rocha escarnitica, parecendo tratar-se de 


ocorrências pobres do ponto de vista econômico não só 
pela irregularidade da mineralização como também pelo seu 
grau de finura. 


Synopsis of articles published in «Técnica» n.º 443 
LII— October 1977 


U, D. C. 621.389:537.311 
H. J. S. ABREU SANTOS 


THE NUMERICAL SIMULATION OF THE MAXWELL-BOLTZ- 
MANN, BOSE-EINSTEIN AND FERMI-DIRAC STATISTICS, 
BASED ON A PSEUDORANDOM NUMBER GENERATOR 


Técnica N.º 443 — LI — 101977, p. 19 a 38 


A comparative study of the obtention of the distribution 
functions corresponding to the Maxwell-Boltzmann, Bose-Eina- 
tein and Fermi-Dirac statistics is done by numerical simula- 
tion on the domain of the small populations of particles. It is 
an application of the Monte-Carlo methods to the simulation 
of a non-determinístic process, that of the random colisions 
of a population of particles, using a pseudo-random number 
generator. Results are presented of the obtention of the equi- 
librium situation for populations of particles of the same 
type, such as Planck's radiation law, as well as for different 
particles, like electrons and phonona. 


U, D. C. 552.1.067.1 


A. F. FERREIRA PINTO 
A. F, BARROS 


OCCURRENCE OF SCHEELITE IN CALCO-SILICATE ROCKS 
FROM VILA FRANCA DAS NAVES AREA 


Técnica N.º 443 — LI|— 104977, p. 39 a 45 


The calc-silicate rocks are concordant with the schists 
which are orientated N 70ºW vertical to sub-vertical. 

These rocks have a granoblastic to porphyroblastic texture 
and they contain plagioclase, diopsidic clinopyroxene, amphi- 
bole and some idocrase and garet. The following accessory 
minerais were also found: sphene, zoisite, clinozoisite, apa- 
tite and opacs. Scheelite was found in small grains in most 
of the rocks studied, 

Some big schist xenoliths occur inside the granite of the 
area and it is thought that the scheelite mineralization and 
the last part of the evolution of the cale-silicate rocks is 
related to the effects of the granite contact metamorphism. 

The scheelite generally occurs in small scattered grains 
in skarnas, so it seems that the economic interest of this 
mineralization will be small. 


= PM SR MME DD Dm 


U. D. C, 378.962.4:550.8[624/628] 


L. AIRES-BARROS 
). F. PINHARANDA 


TEACHING OF GEOLOGY AT |.5.T. TO CIVIL ENGINEERING 
STUDENTS. EVOLUTION AND PROSPECTIVE VISION 


Técnica N.º 443 — LIH|— 104977, p. 1 a 6 


The authors present a summary of the evolution of the 
teaching of Geology to Civil Engineering students at Instituto 
Superior Técnico (Il. S.T.) during the last 20 years. 

After this they present their prospective vision on the 
teaching of Geology and Engineering Geology for that stu- 
dents. The aims of this teaching and the problems regarding 
the need of a balanced knowledge of the fundamentais of 
Geology, their methodologies and their applications to con- 
crete problems of civil engineering are pointed out. The 
main bases considered in order to program the Geology 
course to Civil Engineering students are presented in diagra- 
matic form where scientific and technical fields required for 
the resolution of problem of engineering geology are refered. 

Some comments on examinations are also refered. 


U. D. C. 627.82:627.89 
J. LAGINHA SERAFIM 


RUPTURE OF DAMS DURING 1976. RUPTURE OF TETON 
DAM 


Técnica N.º 44939— LIL— 101977, po. 7a 


This paper concerns the rupture of dams during 1976 and 
continues the work presented by the author in 1960 at the 
LNEC (National Laboratory of Civil Engineering) on the same 
subject. Descriptions of ruptures of relatively small dams 
(above 15m) in North Carolina, Pakistan and Norway are 
presented, Most of the paper concerns the rupture of the 
Teton Dam in the United States in June, 1976. This was the 
highest dam (103m) to ever have failed. The sequence of 
events is presented with photographs. The rupture is attri- 
buted to the very pervious rock, the inefficiency of the grout 
curtains and the existence of a key trench, in which arching 
was produced, possibly allowing hydraulic fracturing. Piping 
and subsequent intense erosion of the soil near the rock 
were certanly major causes. 
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Resumo dos Artigos publicados na «Técnica» n.º 443 
Ano Lil — Outubro 1977 


C. D. U. 624.04 C. D. U. 624.07/012.4])2 | 
ALFREDO DE PINHO MORGADO J. M. MADEIRA COSTA | 
METODO DIRECTO DE DISTRIBUIÇÃO DE MOMENTOS BETÃO NORMAL REFORÇADO COM ARMADURA METÁLICA 
PARA O CALCULO DE ESTRUTURAS SUJEITAS A ESFOR- 
ÇOS HORIZONTAIS Técnica N.c 443 — LI|— 10-1977, p. 51 a 60 


| 
| 
| 
| 


Técnica N.º 443 — LI] — 10-1977, p. 47 a O 


Neste trabalho deduzem-se fórmulas que aplicadas ao 
método de distribuição de momentos em consola reduzem 
este a um ciclo para cada nó da estrutura fazendo desapare- 
cer o inconveniente da lenta convergência, 


= |. 
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Synopsis of articles published in «Técnica» n.º 443 
LI - October 1977 


=. 


U. D.C. eMopozap U. D. C. 62404 
1. M, MADEIRA COSTA ALFREDO DE PINHO MORGADO 
NORMAL CONCRETE REINFORCED WITH STEEL DIRECT MOMENT DISTRIBUTION FOR STRUCTURES WITH 
SIDEWAYS 


| Técnica N.º 443 — Lil — 104977, p. 51 a 60 
| Técnica N.º 443 — LI|— 104977, p. 47 a 50 


The paper works a few formulae to reduce the cantilever 
moment distribution methode to a cycle for each joint of the 
structure with sidesway and overcoming the inconvenience 
of tne slow convergence. 
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1977 


MANUAL DE HIDRÁULICA GERAL 


A. Lencastre, 2.º edição — 1969 


DICIONÁRIO DE UNIDADES E TABELAS 
DE CONVERSÃO 


Vasco Costa e Osvaldo Francês, 1.º edição — 


1959 


NUMEROS E VARIÁVEIS COMPLEXAS 
A. H. Simões de Abreu, 1.º edição 1977. 


450800 


(*) 


250800 


500800 


450800 


75$00 


312550 


(*) LIVROS A PUBLICAR BREVEMENTE 


Desconto de 10% aos assinantes 


GUIA DE ANÁLISE QUÍMICA DAS ÁGUAS 


A. Herculano de Carvalho, 1.º edição — 1961 


CALCULO DE PÓRTICOS 


G. Kani — Trad. da 7.º edição alemã — 1962 


MANUAL DE ESTRUTURAS 
Brazão Farinha — 1.º VOLUME 
Encadernado 


Plastificado 


2º é 3.º Volumes 


BETÃO PREESFORÇADO 


(A cargo dos participantes na Semana de Be- 


tão Preesforçado — |. S. T.— 1975). 


150$00 


100800 


625$00 


587$50 


ei 


(*) 


PEDIDOS À [ERNIES 


Já à Venda 


Números e Variáveis GComplexas 
t EDIÇÃO — 1977 


Eng.º A. H. Simões de Abreu 


Tabelas Técnicas 
8: EDIÇÃO — 1977 


Eng. Vicente Ferreira e Eng. Brazão Farinha 


